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EDITORIAL 


1 — Com a publicação, do primeiro tomo do Traité de Géostatistique Appliqué de G. 
Matheron [7], concretizava-se passo muito importante na Geomatemática — nascia, pelo me- 
nos para os investigadores francófonos, a Geoestatística, que até então tivera incipiente ex- 
pressão em revistas especializadas e relatórios internos de empresas ou serviços públicos. 

Com efeito, quer na versão transitiva, quer no modelo do esquema intrínseco, a Geoes- 
tatística, criada por G. Matheron, mostrava a inconsistência da atitude estatística «ingênua» 
e/ou da posição funcionalista «a todo o custo» que imperavam na Geomatemática. 

Simultaneamente, a Geoestatistica elaborava explicação teórica e, portanto, conferia 
certo conteúdo semântico às evidências experimentais levantadas pelos problemas da avalia- 
ção do teor de blocos e painéis e da evolução da respectiva variância de estimação com a 
dimensão dos mesmos, as quais Krige, Sichel, De Wijs e outros haviam, com rigor exem- 
plar, definitivamente consolidado, uma década atrás. 

Das duas teses de doutoramento de G. Matheron, a primeira Les Variables Rgiona- 
lisées et leur Estimation [8], publicada em 1965, tornou-se pelo rigor e acabamento, fonte 
de inspiração teórica que fez escola, aliás, corporizada no Centre de Morphologie Mathéma- 
tique de Fontainebleau a que se veio juntar o Centre de Géostatistique, ambos ligados à 
Ecole Nationale Supérieure des Mines de Paris. 

Em Fontainebleau, já em 1969, reunia G. Matheron as suas lições no fascículo 1 dos 
Cahiers [11], concatenando aí notável desenvolvimento teórico da krigagem, dentro do modelo 
das funções aleatórias, esclarecendo deste modo, para aquele interpolador linear não condi- 
cional, o seu rigoroso significado e, simultaneamente, explicitando, pela primeira vez, de 
forma correcta, o conceito de deriva (tendência) com o qual a Geofísica e a Geomatemá- 
tica se debatiam, envolvidas em sérias confusões epistemológicas. Não pode, todavia, deixar 
de apontar-se as dificuldades práticas da chamada krigagem universal assim desenvolvida, 
principalmente resultantes do «peso» de cálculo que pode exigir e da impossibilidade, na 
maioria dos casos, de o automatizar [4]. 

Paralelamente, J. Serra, colaborador de G. Matheron, introduz a transformada que tem 
o seu nome e constrói o analisador de texturas que, na actual versão sofisticada, constitua o 
maior avanço da última década no que concerne à morfologia automática. 

A segunda tese de C. Matheron, publicada em 1967 [9] e infelizmente tão pouco 
citada, apesar do recente Random Sets & Integral Geometry [13], influencia fortemente o fas- 
ciculo 3 dos Cahiers [23], onde J. Serra sumariza as suas lições no Centro e que constitua 
marco definitivo para uma tipologia morfológica ainda por acabar. 

Se considerarmos o objecto da Geoestatística — as variáveis regionalizadas, grandezas 
de carácter topo-probabilístico exprimindo casualidade fraca — não é de estranhar que sur- 
jam aplicações, logo na década de 60, em domínios que não o geológico e mineiro, como 
por exemplo, a cartografia automática e submarina, a avaliação e controle da exploração 
de florestas (*), a meteorologia, etc. as quaislevantam problemas difíceis, dispondo-se ape- 
nas, como formalização mais sofisticada, do esquema intrínseco, isto é, das funções alea- 
tórias intrínsecas de ordem 0. 


(1) Segundo G. Matheron dever-se-ia a B. Matem a introdução da terminologia das funções aleató- 
rias no estudo das variáveis regionalizadas (c. f. Materm, B, 1959. Spatial variations. Almaenna Foerlaget, 
Stockoim). 
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Esta limitação é mitigada pela utilização, a partir de 1973, das funções aleatórias intrín- 
secas de ordem k, « que correspondem covariâncias generalizadas, satisfatoriamente interpo- 
láveis, por exemplo, por polinómios de grau não elevado. 

Com efeito, G. Matheron [12] formalizou, sempre no âmbito das funções aleatórias, a 
metodologia apta a tratar variáveis regionalizadas de modo nenhum estacionárias, adaptando 
aos atributos espaciais a conhecida prática usada no tratamento das séries cronológicas 
não estacionárias, a qual consiste em calcular sucessivas diferenças finitas até que surjam 
grandezas quase-estacionárias com as quais se trabalha (*), metodologia esta que vai sendo 
testada nas aplicações levantando, com frequência, e PaRNE operacionais, nomeadamente 
na krigagem universal [14], que é seu dominio 

Por outro lado, problemas práticos derivados do feto da Natureza, no quadro do for- 
malismo das funções aleatórias, apenas fornecer uma realização única daquelas ou da classe 
de funções aleatórias elegida como modelo, levaram a considerar a simulação de variáveis 
regionalizadas nomeadamente no caso de jazidas minerais. 

Assim, em 1974, a tese de doutoramento de A. Journel [5] explicita o trabalho da 
escola de Fontainebleau usando engenhosa metodologia, a qual tem sido cada vez mais apli- 
cada na previsão do valor expectável das respostas das variáveis regionalizadas estacioná- 
rias a várias classes de operadores, embora a sua utilização aos casos nitidamente não 
estacionários só tenha sido esboçada muito recentemente [14]. 

O impasse da realização única, incluso na conceptualização que recorre às funções 
aleatórias, obteve, deste modo, uma saída pragmática rigorosa no seu próprio âmbito. 

A estrutura revelada até ao momento de ordem dois e pelo histograma, na única reali- 
zação, constitui o condicionamento da simulação, a qual tem, portanto, o valor que tiver 
a inferência possível a partir de tal realização. Porém, a não existência de testes de signifi- 
cância para a prática da inferência, ou a impossibilidade de usar aqueles de que se dispõe, 
leva, por vezes, os estatísticos ortodoxos a olhar suspeitosamente aquelas simulações e, dai, 
que, aqui também e cada vez mais, se fale da robustez das estatísticas inferidas e, em deses- 
pero de causa, se multipliquem as simulações (g raças ao recurso ao computador ) para mostrar 
que as esperanças sobre o conjunto daquelas é coerente com as médias espaciais obtidas, 
por exemplo, por krigagem, explicitando-se, deste modo, por via «experimental», à ergodici- 
dade implícita no próprio modelo. 

Entretanto, multiplicam-se por toda a parte as aplicações e, por força destas, surge a 
Geoestatística mineira onde avulta a revisão do conceito de parametrização, com a sua base 
teórica alicerçada em certo tipo de convexo que é, mais uma vez, G. Matheron que estru- 
tura em 1975. 

E, como não podia deixar de ser, «software» é intensamente desenvolvido para apl 
cação a casos concretos. 

Para lá do problema da inferência, a Geoestatística mineira enfrentou duas dificul- 
dades igualmente importantes, depois de haver posicionado com rigor o significado do suporte. 

Um deles é o da impossibilidade de conhecer, excepto em alguns casos favoráveis, o 
histograma dos valores reais da varidvel regionalizada, uma vez determinado o histograma 
dos interpoladores de krigagem (que é um estimador linear não condicional e não envie- 
zado), embora se saiba encontrar a variância de dispersão das grandezas reais. 

O outro consiste em prever qual será o histograma dos futuros estimadores em fase 
adiantada do reconhecimento mineiro ou mesmo durante a exploração da jazida, havendo 
disponível mais informação, o que é de importância evidente, porquanto será à custa dessa 
informação adicional que se tomarão, na maioria dos casos, as decisões ao nível da exploração. 

Este segundo problema conduziu ao aparecimento da Geoestatística não linear onde é 
relevante o conceito de função de transferência introduzido por G. Matheron, à custa do qual 
é possível realizar a krigagem disjuntiva, isto é, construir um estimador condicional fraco 
dos teores dos blocos sobre os quais se irão tomar as opções de exploração e tratamento. 

Notar, todavia, que a krigagem disjuntiva fornece apenas os histogramas expectáveis 
tendo por suporte pequenos blocos. 

Por outro lado, convém referir que, nas aplicações, o variograma se tem tornado o 
instrumento universal da Geoestatística linear. Os seus utilizadores empregam-no para a inter- 
pretação estrutural das variáveis regionalizadas que estudam. 


E 


(1) De referir que o conceito de variograma, tão usado em Geoestatística, também surge, ao que parece 
pela primeira vez, no estudo das séries temporais (c.f. Jowett, G. H, 1954. Least squares regression analy- 
gis for trend-reduced time series. J. Roy. Stat. Soc., Series B 17, 91-104). 
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Esta tendência é necessário que se assinale, porque traduz, por vezes, a tentativa de 
esboçar tipologias baseadas em relações estruturais. 

E óbvio que o uso sistemático do variograma como revelador da estrutura levanta 
problemas intrincados, principalmente no que respeita à escolha do seu modelo teórico. 

Com efeito, esta escolha está quase sempre ligada a uma interpretação da génese e/ou 
metalo-génese da mineralização, etc., o que confere à análise variográfica, isto é, à aplicação 
prática do variograma, aspectos determinísticos exteriores, evidentemente, à metodologia das 
funções aleatórias, que tem sido o suporte da Geoestatística. 

Contudo, esta dificuldade não surge quando o variograma é apenas usado com vista à 
realização da krigagem, porque neste caso o estimador correspondente é relativamente 
robusto, isto é, não é sensivelmente afectado pelo detalhe do modelo adoptado para o vario- 
grama. 

Ainda no aspecto da análise estrutural, aparecem dificuldades relacionadas com a 
forma de construir o próprio variograma experimental. Na verdade, a Geoestatística parte 
da hipótese de se dispor de amostragem regular, o que nem sempre acontece, como é sabido. 
Dai que tenha aparecido a prática de utilizar variogramas regularizados [3] e mesmo se 
comecem aq usar médias ponderadas para efectuar o seu cálculo. 

Além disso, quando existe deriva, a definição do variograma das flutuações é passiva 
de várias interpretações criadas com a finalidade, talvez falaz, de encontrar «objectiva- 
mente», duas coisas simultaneamente: a deriva .e o variograma residual. 

Aqui também, o modelo do variograma é crucial. Todavia o problema deixou de ter 
a relevância tradicional que auferia na Geomatemática, graças à utilização das covariâncias 
generalizadas, decorrentes das funções aleatórias intrínsecas de ordem kh. 

2 — Entre os portugueses, a Geoestatística é aplicada pela primeira vez por J. Rogado 
em 1964 [17] no cálculo da reserva de jazidas tabulares de minérios de ferro de Cassinga 
(Angola ). 

Como seria de esperar, a utilização da metodologia da Geoestatística levantou celeuma. 
Todavia, este facto não impediu que se procedesse, de forma sistemática, à avaliação das 
reservas dos jazigos eluviais de minérios de ferro da área referida, utilizando a Geoestatis- 
tica, tarefas que terminou em fins de 1966 [18]. 

Porém, o engrossar da celeuma levou à deslocação a Angola de G. Matheron, que con- 
firmou a justeza do cálculo usado [1]. 

De passagem por Lisboa, G. Matheron efectuou uma conferência no IT. 8. T., da qual 
resultou a publicação de um artigo na revista «Técnica» [10]. 

Apresentava-se, assim, no foro universitário, a Geoestatística já praticada, em Portu- 
gal, por docentes do I. 8. T. 

Há todavia que sublinhar ter sido a necessidade objectiva de avaliar correctamente 
importante reserva mineral que suportava um projecto concreto, que tornou possível a inser- 
ção eficaz deste ramo de conhecimento na Universidade. 

Assim, continuando as aplicações da Geoestatística, tornou-se evidente a necessidade de 
desenvolver a formação, neste domínio, quer de docentes e investigadores universitários, 
quer de engenheiros colaborando no processo industrial. Dai que, já em 1969, dois docentes 
do I. 8. T. tenham frequentado cursos em Fontainebleau e que, nos anos seguintes, mais 
docentes e engenheiros continuassem com esta prática. 

Durante 1970, realiza-se, no I. 8. T., o primeiro curso livre de Geoestatística leccionado 
pelo Eng.º L. Cortez; em 1971, publica-se na «Técnica» o primeiro artigo em lingua portu- 
guesa reportando uma aplicação prática neste domínio [16] ('). 

Em 1971 é introduzida no «curriculum» do curso de Minas do TI. S. T. a cadeira de 
Geoestatística, que em 1975 se desdobra em duas outras e em 197h a A. E. I. 8. T. edita o 
primeiro curso de Geoestatística [15]. 

No princípio da década de 70 funciona no Núcleo de Tratamento de Minérios do I. A. C., 
sedudo no I. 8. T., um seminário de Geoestatística, cujos colaboradores mantêm contactos 
frequentes com o Centro de Fontainebleau, contactos esses que têm continuado dentro do 
âmbito do C. V. R. M. U. T. L., herdeiro da investigação do referido Núcleo. 

Simultaneamente, continuam as aplicações da Geoestatística graças à colaboração esta- 
belecida entre o TI. 8. T. e a Companhia Mineira do Lobito, com minas em Angola ea J. E. N. 
com minas em Portugal. 

Nasce, assim, a escola de Geoestatística de Lisboa que, aliás, é a primeira resultante 
de Fontainebleau fora da França. 


O Co 


(1) De notar que em 1969 é editado texto sobre axiomática da Geoestatística [6]. 
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A escola de Lisboa desenvolve principalmente métodos de Geoestatística mineira. 

Neste domínio avultam duas linhas de desenvolvimento, na fase inicial de trabalhos 
do grupo de Lisboa, a saber: 

(i) — A análise vario-grupal é posta em forma prática e as suas virtualidades 
testadas em minérios de ferro e de uránio. 

Trata-se de método de análise estrutural que utiliza uma janela para escrutini- 
zar a estrutura das mineralizações recorrendo, depois, à lógica vaga para classificar 
os resultados obtidos. Seguidamente, introduz uma noção de proximidade entre as 
estruturas assim reveladas, a qual permite, realizar tipologia e orientar problemas 
sequenciais de pesquisa mineira. 

(ui) — A parametrização das jazidas e o seu relacionamento com as políticas 
de rejeição são vistos na escola de Lisboa principalmente em termos de variáveis 
mineralúrgicas, sugerindo o abandono do teor tal-qual como parâmetro básico de 
controle do projecto mineiro. 

Neste domínio, desenvolve o grupo de Lisboa metodologia peculiar que permite, 
à custa da análise convexa, estabelecer políticas optimais de exploração e tratamento 
de jazidas minerais. 

Partindo duma ideia inicialmente desenvolvida por R. Tillé [24], constroe-se, no 
domínio do discreto, o convexo superior que relaciona a recuperação em metal com 
o rendimento em peso de concentrado |20, 21, 22], o qual tem como caso particular 
o convexo superior ligando o metal tal-qual com o peso do minério tal-qual extraído, 
que havia sido elaborado por Fontainebleau. 

Assim, a metodologia de parametrização desenvolvida em Lisboa toma em con- 
sideração aspectos do processo mineralúrgico que aquele outro convexo não pode cobrir. 

Como seria natural, a Geoestatística mineira praticada pelo grupo de Lisboa influen- 
cia fortemente o Planeamento Mineiro que surge paralelamente como opção dentro do Curso 
de Engenharia de Minas do TI. 8. T. 

É neste domínio que a escola de Lisboa tem ultimamente mais insistido, elaborando 
metodologia que permite a programação de desmontes, quer a céu-aberto, quer subterrâneos, 
a qual se traduz no constituição de poderosos «packages» de «software» que têm sido utili- 
zados em aplicações concretas importantes, em jazidas nacionais e estrangeiras. 

Outra preocupação da escola de Lisboa, trata-se da taxonomia automática (análise gru- 
pal, análise factorial, análise das correspondências, etc.). 

Aliás, segue-se aqui a mesma tendência que se manifesta um pouco por toda a parte 
entre os cultores da Geociência. 

Também facto que não deve ser esquecido resultante da actividade da escola de Lisboa 
é a recente criação, no Brasil, de um núcleo de Geoestatística, encabeçado por antigo aluno 
do TI. 8. T., o qual surgiu após regência de um curso de Geoestatistica mineira na Universi- 
dade de S. Paulo, em 1976, por J. Rogado, o qual tem realizado, naquele país, trabalho de 
aplicação muito relevante com significativos contactos com Lisboa. 

3-— A teoria da amostragem de materiais a granel de P. Gy permitiu, pela primeira 
vez neste domínio, satisfazer simultaneamente exigências teóricas mínimas e dar resposta 
eficaz em aplicações frequentemente complexas. 

Entre nós, a metodologia de P. Gy tem vindo a ser utilizado e desenvolvida desde a 
década de 60, a nível laboratorial e industrial, nomeadamente no Departamento de Minera- 
lurgia e Planeamento Mineiro do I. 8. T., no reconhecimento de jazidas tabulares de miné- 
rios de ferro de Cassinga (Angola), no controle da sua exploração e tratamento e ainda na 
concepção da instalação automática de amostragem do Porto de Moçâmedes. 

O desenvolvimento das suas aplicações e inserção no «curriculum» duma cadeira do 
Curso de Minas do T. S. T. conduziu à necessidade de implementar a formação, neste domi- 
nio de investigadores universitários e de docentes colaborando na sua aplicação industrial. 
Dai que, logo em 1971, dois docentes do I. S. T. tenham frequentado cursos em Cannes, regi- 
dos por P. Gy. 

Na década de 70 continuam a realizar-se aplicações da teoria, graças à colaboração dada 
à Companhia Mineira do Lobito e, ultimamente, à Comissão para o Lançamento do Programa 
do Aproveitamento Integrado das Pirites Alentejanas e Pirites Alentejanas E. P. 

Assim, à teoria e a prática da amostragem no âmbito do modelo gizado por P. Gy, 
que é, aliás, uma aplicação sofisticada da Estatística, passou também a ser preocupação da 
Escola de Lisboa, pois o seu objecto diz respeito a variáveis regionalizadas ou a séries 
temporais. 
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4 — Em seguimento, apresentam-se a maioria das intervenções e as conferências do 
Beminário pela ordem por que foram dadas e espera-se que a sua publicação possa contri- 
buir para a divulgação da Geomatemática nos países de língua portuguesa. 

Faz-se notar que este Seminário teve intencionalmente carácter introdutório, pelo que, 
embora as contribuições nele apresentadas se situem em domínios avançados da Geomatemá- 
tica, não houve «a preocupação de apresentar resultados que não fossem já conheculos ou 
que não estivessem publicados. 

Por outro lado, também intencionalmente, incluiu-se uma conferência sobre modelos 
mineralúrgicos como sugestão para a construção da ponte necessária à implantação correcta 
das variveis mineralúrgicas nos conceitos topo-probalísticos da Geoestatística. 

5— A Comissão Organizadora do Seminário e Conferências de Geomatemática agra- 
dece calorosamente a contribuição dos autores nacionais e estrangeiros que vieram enrique- 
cer, quer o Seminário, quer as Conferências. 

Também a Comissão Organizadora deseja expressar o seu agradecimento aos participan- 
tes nacionais e estrangeiros, nomeadamente aos docentes da Faculdade de Ciências da U. L. 
e da Faculdade de Ciências e Tecnologia da U. C. e aos colaboradores da Direcção-Geral de 
Minas, do Serviço de Fomento Mineiro, da Junta de Energia Nuclear, do C. P. R. M. e. P. 
T. do Brasil e das empresas Paulo Abib Engenharia S/A do Brasil, Ferrominas E. P., Pirites 
Alentejanas E. P. e Mineira de Santiago, 8. A. R. L. cuja presença interessada e actuante foi 
a maior compensação para os esforços realizados na organização deste Seminário e Confe- 
rências. 

O Laboratório Nacional de Engenharia Civil e a Sociedade Portuguesa Time-sharing 
8. A. R. L. são merecedoras do nosso reconhecimento pelas facilidades de cálculo automático 
postas à disposição do Seminário. 

Uma palavra de apreço e agradecimento deseja ainda a Comissão Organizadora ende- 
reçar aos alunos do 5.º ano do Curso de Engenharia de Minas de então e aos colaboradores 
do C.V. R.M.U. L. e do Departamento de Mineralurgia e Planeamento Mineiro do I. 8. T. e 
bem assim a todos os que actuaram na dactilografia e pre-iimpressão dos trabalhos, sem os 
quais teria sido impossível levar ávante o Seminário e Conferências. 

Em especial, a Comissão Organizadora deseja manifestar a sua gratidão ao Eng.º Fer- 
nando Soares Carneiro, Director-Geral de Minas pelo apoio concedido e, bem assim, por ter 
participado nos trabalhos do Seminário na sessão de abertura, 

Particular reconhecimento quere ainda a Comissão Organizadora manifestar à Presidên- 
cia do I. N. I. €C. pelos fundos concedidos, quer para a organização do Seminário, quer para 
ajuda da edição que agora se efectua. 


J. QUINTINO ROGADO 
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Introdução às variáveis regionalizadas 


RESUMO 


Introduz-se a noção de Variável Regionalizada 
no seu duplo aspecto aleatório e estrutural. Dão-se 
exemplos de Variáveis Regionalizadas nos dominios 
geológico e mineiro e sugere-se, em Referências Bi- 
bliográficas existentes no CYVRMUL, outras aplica- 
ções. Interpretam-se as Variáveis Regionalizadas 
como realização de uma Função Aleatória e introdu- 
zem-se hipóteses restritivas. Define-se o variograma 
como instrumento básico da Geoestatística e estulam- 
-se algumas das suas propriedades como revelador de 
estruturas. 

Indica-se o método de cálculo do variograma 
experimental e listam-se alguns modelos teóricos 
ajustáveis. 
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1 — GENERALIDADES 


Nas Ciências da Terra surge um tipo de variá- 
veis distribuidas no espaço e que apresentam uma 
certa estrutura. 

Se for x a coordenada (*) do ponto do espaço 


(*) Em tudo o que ge segue, x é a coordenada a 1, 2 
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ABSTRACT 


The theory of Regionalized Variables is introdu- 
ced by its structural and random. features. Examples 
of the variables in the geological and mining fields 
are presented and other applications are suggested. 
The Regionalized Variables are studied as a realiza- 
tion of Random Functions and work hypothesis are 
introduced. The variogram is definned as the basic 
tool of Geostatistics und some of its applications as 
structure depictor are discussed. 

The method for experimental variogram calcula- 
tion is described and some theoretical models are 
udjusted. 


a que se refere a variável, Y(x) designa o valor que 
a variável toma nesse ponto. 

Por exemplo, o teor ao longo de uma sondagem 
é uma Variável Regionalizada (V, R.). 

Estas variáveis apresentam-se segundo um duplo 
aspecto contraditório: têm uma característica aleató- 
ria que se manifesta na irregularidade e variação 
imprevisível de um ponto para outro e uma caracte- 
ristica estruturada que reflecte as ligações existentes 
entre os pontos do espaço onde o fenómeno minerali- 
zado ocorreu. 

As variáveis regionalizadas não podem ser trata- 
das pela estatística clássica (Bernoulliana) visto que 
esta implica que: 


1. As provas possam ser indefinidamente repe- 
tidas 


2. As provas são independentes (o resultado de 
uma prova não pode ser influenciado pelo 
das que a precederam). 


Ora estas condições não são respeitadas nas va- 
riáveis tratadas nas Ciências da Terra: tirando uma 
amostra num ponto de coordenada x, o teor de tal 


ou 3 dimensões. 
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amostra é um valor único, fisicamente determinado, 
e não é possivel voltar a repetir a experiência. No 
entanto, se tirarmos outra amostra num ponto infini- 
tamente vizinho, poder-se-ia dizer que a condição 1. 
era respeitada numa certa vizinhança, mas a condi- 
cão 2. deixava de o ser, visto que, se a primeira 
amostra fosse tirada numa zona rica, havia tendên- 
cia, em média, para que a segunda tivesse também 
um teor elevado. 

Por outro lado, um aspecto fundamental das 
V. R. que é a continuidade do fenómeno mineralizado, 
isto é, o modo mais ou menos regular como estas 
evoluem no espaço, não pode ser tomada em conta 
pela estatística clássica. 

De facto, consideremos o seguinte exemplo de 
amostras de solo retiradas ao longo de 2 linhas a 
intervalos regulares e suponhamos que a variável 
que nos interessava era a percentagem de água 
nessas amostras. Os resultados são indicados na ta- 
bela seguinte: 


5 10 15 20 25 20 15 10 5 
10 2 15 10 20 5 15 5 20 


linha A 
linha EB 


Do ponto de vista da estatística clássica, as 
amostras tiradas na linha A não diferem das tiradas 
na linha B — têm a mesma média, variância, momen- 
tos de ordem superior, histograma, etc. No entanto, 
do ponto de vista da distribuição do teor da água 
no solo é intuitivo que as duas sequências são pro- 
fundamente diferentes e portanto há um carácter nas 
V. R. que tem de ser tomado em conta através de 
outra ferramenta que não a estatística clássica. 

Algumas características qualitativas das variá- 
veis regionalizadas são comentadas seguidamente: 

a) Localização — a V, R. estã definida numa 
determinada região do espaço à qual se chama cam- 
po geométrico que pode ser uma formação geológica, 
por exemplo, um jazigo em que as variáveis têm um 
conjunto de ligações entre elas dado pela própria 
génese (fenómeno estruturante). Mas na prática a 
V. R. não é uma função pontual — é o valor médio 
numa amostra de dimensões finitas, perfeitamente 
determinadas. O domínio onde a V. R. é medida de- 
signa-se por suporte. Para um teor, por exemplo, o 
suporte é o volume da amostra onde este foi anali- 
sado. Se se muda o suporte, obtém-se outra regiona- 
lização diferente da primeira até do ponto de vista 
da estatística clássica. O histograma das carottes 
de sondagem apresentará sempre dispersão maior 
(maior variância) do que o histograma de blocos de 
algumas toneladas (Fig. 1). 


pr 


b) Continuidade —- é a variação espacial da V. 
KR. Apesar da complexidade das flutuações verificadas 
existe em geral continuidade em média. Esta defi- 
ne-se por 

lim E(Y(x)-Y(x ))—o 
XX, 

Há casos em que tal continuidade não se verifica 
e então o limite C £o para o qual tende a diferença 
dos quadrados da variável quando os pontos se apro- 
ximam indefinidamente é o chamado efeito de pepita. 
Dc facto, na fronteira de uma pepita existe descon- 
tinuidade: 


Pepita Esteril 


Xx— X5 


] 
Lim [Y00)-Y()]"= 0,%0 
O 

c) — Anisotropias — em certos casos, ligados a 
fenómenos geológicos particulares, a continuidade 
numa determinada direcção pode ser diferente da 
que se verifica noutra direcção, Ao longo de uma 
mesma camada sedimentar é de esperar uma varia- 
ção muito mais regular do que na direcção perpen- 
dicular. 

d) Fenómenos de transição — No interior do 
campo geométrico da variável, podem surgir estrutu- 
ras particulares como sejam conjuntos lenticulares 
em que existe continuidade no interior de cada uni- 
dade e variação brusca de unidade para unidade 
(Pig. 2). 

Este conjunto complexo de descontinuidades rea- 
liza uma partição do espaço mineralizado em compar- 
timentos mais ou menos independentes (*). É de es- 


Fig. £ 


(*) Corresponde, por exemplo, a micro-bacias de sedimentação onde as condições de deposição evoluíram de um 


modo mais ou menos autónomo. 


90 


TÉCNICA 451/452 


—em... ie im cm e e um e TT a TS —— 


- Mm) a 
4 A 
//H 


“Postes de cimento 


armado Cavan 


Qualidade que 
desafia o tempo 


Av. Visconde Valmor, 76-1.º 


- Tel. 766014 (7 linhas! Lisboa-l 


“Ventilação 
de Edifícios 
Ventilação Térmica FÉRIA 
(tipo shunt) 
ECONOMIA E EFICÁCIA! 


Com aprovação do 
Laboratório Nacional 
de Engenharia Civil 
Câmara Municipal de Lisboa 
Batalhão de Sapadores 

E Bombeiros 
Diracção-Geral de Saúde 


Consulte o Ar 
das Editicaçõe 


é Respirar! 


FÉRIA SERIC 


ua É Sebastião da 
dreira 10 18 | 


Tels 54 938º 5769771 | 


Consulte O nosso 
Gabinete Tecnico 


ESPECIALISTA DESDE 1947 


, sopecate 


ESTUDOS 
GEOTÉCNICOS 


FUNDAÇÕES 
ESPECIAIS 


RUA DO ARSENAL, 146-2º 
TELFS.320208:360437:364010 
LISBOA 


SONDAGENS 
Geológicas 

Mineiras 

Geotérmicas 

De Captação de águas 


FUNDAÇÕES 
TRABALHOS MINEIROS 


Abertura de poços e galerias 
Injecções 
CONGELAÇÃO DE SOLOS 


Para trabalhos mineiros 
e de construção civil 


R. Cruz da Carreira, 6, 3.º 
Tel. 52797 
Lisboa 1 


fora 


TECNICA ll 


Mais do que 


m milhão de HP 


na Europa 


mas, nem um HP a mais! 


A potência de 
todas as tur- 
binas monta- 
das nas má- 
quinas +6F+ 


de decapagem a jacto de grenalha ultra- 


passou o primeiro milhão de HP, o que é 
muito, mas, não é demais a energia econo- 
mizada pela força centrífuga desse 
mesmo jacto. 

No entanto, o proveito obtido, em relação 
às superfícies racionalmente limpas de 


metal, plástico 
ou borracha 
contam mais 

de que os va- 

e sá SN lores em HP, 

Seria bom se fosse possível medir, pelo 
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perar que o sistema de relações entre a V. R. no 
interior de cada lenticula seja mais «forte» do que 
quando se passa de uma para outra. Esta estrutura 
que conduz a que, a partir de uma certa distância, 
haja independência entre as V. R., é designada por 
Fenómeno de Transição. 


+ 
* + 


A teoria das Variáveis Regionalizadas foi desen- 
volvida a partir dos trabalhos de Matheron [1], [21]. 

A aplicação desta teoria a problemas mineiros, 
em especial cálculo de reservas e planeamento da 
exploração, designa-se por Geoestatística. 

Um resumo do formalismo e campos de aplica- 
cão da Geoestatística encontram-se em [8], [4], [5], 
[6] que constituem guias práticos de utilidade para 
um primeiro aproach dos problemas levantados pela 
aplicação da teoria das Variáveis Regionalizadas no 
domínio mineiro. 

Considerações gerais sobre a utilidade da Geoes- 
tatística do ponto de vista da empresa mineira 
(utilização como instrumento de decisão económica) 
encontram-se brevemente explanadas em [7], [8], 
[9] e [101. 

Como exemplo de Variáveis Regionalizadas, te- 
remos, como primeiro caso clássico, o teor num de- 
terminado elemento numa amostra localizada no es- 
paço de coordenada x. De facto, o teor apresenta 
uma evolução errática de ponto para ponto do espaço, 
mas será de esperar que teores em pontos vizinhos 
tenham mais relação entre si do que teores em pon- 
tos muito afastados. Ao aspecto aleatório, sobrepõe-se 
a estrutura do fenômeno mineralizado que dá uma 
certa unidade ao conjunto de variáveis medidas no 
campo geométrico que é a jazida. 

A possança mineralizada é também uma V. R. 

Variáveis de tratamento como rendimentos são 
também V. R. pois a aptidão ao tratamento depende 
de características intrínsecas do minério que se de- 
senvolvem no espaço de uma certa maneira, depen- 
dente do fenómeno mineralizado. 

Além destas variáveis «mineiras» (cuja aplica- 
ção no cálculo de reservas é imediata) existem mui- 
tas outras que podem ser tratadas como V. R. 

Por exemplo, o estado de fracturação de um 
mármore, definido como a superficie total de fractu- 
ras por unidade de volume é objecto de um trata- 
mento geoestatístico em [111]. 

Na prospecção geoquimica, a dimensão e varia- 
bilidade das amostras colhidas para estudar a dis- 
tribuição de determinado elemento no espaço (homo- 
geneidade, dispersão, anisotropias, técnicas de amos- 
tragem e «mapping») são aplicações típicas da teoria 
das V. R. [12]. 

Em hidrogeologia e pluviometria, a altura piezo- 
métrica num determinado ponto de um aquífero ou 
a quantidade de precipitação que cai num determi- 
nado ponto de medida podem ser tomadas também 
como V. R. [18], [14]. 

Em geofísica, em particular na interpretação de 
dados gravimétricos [15], mas também em meteoro- 
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logia, em que se pretende tomar em conta a relação 
entre a pressão e a velocidade dos ventos [14], a 
teoria das V. R. tem-se revelado como instrumento 
poderoso de interpretação e previsão. 

A determinação de cartas dos fundos submarinos 
é também um caso em que a cota medida num deter- 
minado ponto pode ser tomada como V, R. [14]. 

Em problemas de controle da poluição [16] e 
outros um pouco fora do âmbito deste trabalho (den- 
sidade de árvores numa floresta, distribuição geográ- 
fica da população e actividades econômicas) também 
a teoria da V. R. pode ter aplicação útil. 

EM RESUMO, qualquer variável dependente do 
espaço, em que, para além de um. carácter aleatório, 
existe um sistema de relações entre os valores toma- 


dos em pontos distanciados de um vector q (e esse 
sistema de relações não pode ser interpretado de um 
modo simples a partir de uma função que descreva 
uma realidade «física») pode ser tratada como Farid- 
vel Regionalizada e sofrer uma análise segundo o 
formalismo desenvolvido pela Geoestatística. 


2-— AS VARIÁVEIS REGIONALIZADAS COMO REALI- 
ZAÇÃO DE UMA FUNÇÃO ALEATÓRIA 


Para interpretar o duplo aspecto aleatório e es- 
trutural das V. R. utiliza-se a noção de Função Alea- 
teória (F. A.). 

Uma F. A. é uma variável aleatória a uma infi- 
nidade de componentes. 


Nejata %,, Xi, cesso , X, as coordenadas dos pontos 
do espaço onde se realizam medidas das V, R. Os va- 
lores T(x,)), T(X,), css » Y(x,) — valores que 


a V. R. toma nesses pontos — podem considerar-se 
como uma realização particular de uma F. A. am 
componentes, isto é, como um. dos resultados possíveis 
de uma tiragem à sorte da F. A. Se fosse possível 
repetir a experiência, obter-se-ia outras realizações 
da mesma F. A. cuja lei de distribuição de probabi- 
lidade a n componentes seria f[T(x)), TX), e... i 
Y(x,)]- 

Como a partir de uma única realização de uma 
variável aleatória ordinária (um número tirado ao 
acaso) não é possível inferir a lei de distribuição 
dessa variável, também aqui a lei de ditribuição de 
probabilidade a mn componentes não se pode determi- 
nar, 

Exige-se assim a introdução de hipóteses restri- 
tivas. 


4 — HIPÓTESES RESTRITIVAS 


Para ultrapassar a dificuldade acima enunciada 
podem formular-se algumas hipóteses restritivas, 

a) Estacionaridade — A F, A, é estacionária se 
a sua lei de distribuição de probabilidade é invariante 
por translacção, isto é a F. A. 

[EAR EM rioãs » Y(x,)] tem a mesma lei 
de distribuição de probabilidade do que 


[Ex th), TOrtTh), css S T(x, +h)] 
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neste caso E(Y(x)j=m(x)=E(Y(x+h)=m=ce'!* 
portanto a média é invariante por translacção e pode 
ser tomada como nula, . 

b) Estacionaridade de 2.º ordem. 

Se a covariância 


K(h) =E(Y(x+h). Y(x)j-m(x)? 


não depende de x (isto é, do ponto de apoio) e só 
da distância que liga os dois pontos entre os quais 
é calculada, diz-se que há estacionaridade de 2.º 
ordem. 

Esta hipótese implica a existência da variância 


K(o) =E(Y(x)2)-m(x)? 


finita, o que nem sempde se verifica (*). 

c) Hipótese intrínseca. 

Consideremos os incrementos [Y(x+h)-Y(x)] 
e admitamos que estes são estacionários, isto é, não 
dependem de x. 

Então 
E([Y(x+h)-Y(x)]j=m(h) deriva linear (que se 
pode supor nula) 
E([Y(x+h)-Y(x)]2)=2 y(h) variograma 

No caso da hipótese intrínseca, o variograma 
y(h) é independente do ponto de apoio e portanto é 
suficiente para definir a estrutura da V. R. no espaço 
em que a estacionaridade se verifique. 


Estas hipóteses a), b) e c) são sucessivamente 
mais fracas e devem ser testadas para investigar se 
num determinado campo se verificam. E esta preci- 
samente a dificuldade inicial do técnico de geoesta- 
tística aplicada — determinação de um variograma 
único para todo o espaço que pretende tratar (ou a 
partição do espaço em zonas em que a estacionari- 
dade, tal como se definiu para a hipótese intrínseca, 
seja válida). 

De facto, para que o fenómeno geológico admita 
a hipótese intrínseca, é necessário que o sistema de 
correlações entre pontos do espaço só dependa da 
distância e não da zona de trabalho (medida pela 
coordenada x). Isto implica que o fenómeno geoló- 
gico tenha agido de um modo «aproximadamente» 
constante para a zona onde se pretende estimar o 
variograma e que a escala da rede de amostragem 
«acerte» com a escala em que o fenómeno revela a 
estacionaridade (*). 

Mas há fenómenos que, seja qual for a escala 
de trabalho, nunca apresentam estacionaridade. Um 
fenómeno governado pela gravidade (o teor aumenta 


sempre em profundidade, por exemplo) não pode evi- 
dentemente sofrer o tratamento da geoestatística es- 
tacionária. Para tais casos, a Krigagem Universal 
[18] é a metodologia a utilizar. Vamos no entanto 
cingir-noz, meste trabalho, ao caso em que a estacio- 
naridade, mesmo local, pode ser detectada. 


4 — VARIOGRAMA 
4.1 — Introdução 


O variograma é o instrumento básico da geoes- 
tatística. É a partir desta função que todas as esti- 
mações são efectuadas e afectadas de um erro. Por 
outro lado, é um poderoso revelador de estruturas. 
É neste último aspecto que o variograma vai ser 
considerado neste trabalho. 

Num determinado campo geométrico, considere- 
mos o valor que a variável regionalizada toma em 


dois pontos x e x + onde h é um vector. 

Na hipótese intrínseca, o variograma não depen- 
de do ponto de apoio e pode ser calculado para qual- 
quer região do espaço. 


É portanto o valor médio do quadrado das dife- 
renças entre todos os pares de pontos existentes no 


campo geométrico distanciados de n Ê 

Na prática, é necessário que o número de pontos 
amostrados seja suficiente para que a média tenha 
significado estatístico. 


4.2 — Fenómenos de transição 


O variograma, nos fenómenos de transição, apre- 
senta-se com o aspecto da figura seguinte (vario- 
grama do teor em Cu num jazigo sedimentar): 


10 20 30 —h 


Variograma dos teores em Cu (segundo [19]) 


(*) Um estudo experimental realizado por Krige [9], [17] verificou que a variância de uma amostra num campo 
sucessivamente crescente (até atingir todo o jazigo do Rand), aumentava indefinidamente segundo gz(v/v)=a log (v/v) e 
portanto não há garantia de que exista variância à priori finita. Sendo py a dimensão da amostra e v a dimensão do campo 


onde a variância g2(v/v) é calculada, verifica-se que lim 


Vs 


o2(v/v) = 00. 


1&s 


(*) É neste aspecto que a informação geológica é muito útil na formulação dos modelos geoestatísticos. A interdepen- 
dência entre as observações de campo e o tratamento matemático deve ser constante, 
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Na maioria dos casos práticos, há um ramo 
crescente do variograma para valores pequenos de 


Nes 
he a partir de uma certa distância, este tende para 
um patamar. Este facto é interpretado teoricamente 
a partir da estacionaridade de 2.º ordem: 


)=2 E À [esmo 
) = ZÉ [rea + h)]º ! -— E [r(x +). Y(x)| ci 


| I 
TE LO] 
em virtude da estacionaridade de 2.º ordem, 
E([Y (x + n)]2) = E(UY(x)]) = K(o) 


y(h)=K(o) - K(h) 


Nos fenómenos de transição, a partir de uma 
certa distância a, dita amplitude, a covariância K(h) 
anula-se e o variograma reduz-se ao patamar € que 
coincide portanto com a variância à priori. 

Outro parâmetro característico do variograma é 
o efeito de pepita C que corresponde à ordenada na 
origem. 

E o caso que se apresenta na figura seguinte, 
correspondente às espessuras mineralizadas de um 
jazigo sedimentar de Ferro: 


O efeito de pepita pode ser interpretado como 
o «testemunho» de micro-regionalizações existentes 
a escalas inferiores à de amostragem. De facto, o 
primeiro ponto a que temos acesso experimental é 
y(1) e não se sabe o andamento do variograma entre 
os pontos O e 1; apenas se pode extrapolar linear- 
mente para o eixo das ordenadas. 

Suponhamos o caso ilustrado na Fig. 3. 


yth) | 


Fig. 3 
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O variograma A foi construído a partir de 
amostras cuja distância mínima era 1. Se for aper- 
tada a malha de amostragem, poder-se-ia revelar o 
variograma B que corresponde à sobreposição de uma 
transição a pequena escala (com o seu patamar) 
com a transição à escala da 1." amostragem. Além 
deste aspecto («testemunho» de microregionaliza- 
ção), o efeito de pepita reflecte ainda os erros de 
amostragem. De facto, se fosse possível realizar a 
colheita de outra amostra no mesmo ponto x, essa 
amostra daria um valor para a V. R. Y' diferente 
do inícial devido aos erros de amostragem, analíticos, 
etc. Portanto, duas amostras infinitamente próximas 
podem contribuir para o termo [Y(x+h)-Y(x)I> 
apenas pelos erros de amostragem. 

Outro aspecto importante da mineralização reve- 
lado pelo variograma é a presença de estruturas 
embricadas e corresponde ao caso ilustrado na Fig. 5 
para o variograma global da acumulação (teor x 
espessura) de um jazigo de Au, 

Este caso corresponde ao facto de a malha de 
reconhecimento utilizada ter «surpreendido» simulta- 
neamente duas estruturas de transição embricadas 
segundo o esquema da Fig. 4: 


Fo =] Das... o = = 


Sd mM 12 193 ntoZm 
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Fig. 5 


O problema da escala encontra-se aqui de novo: 
para pequenas distâncias, o variograma reflecte a 
transição a essa escala (no interior das lentículas) e 
para distâncias superiores à dimensão média das 
lentículas, o variograma reflecte o modo como estas 
se estruturam entre si. 


Em resumo, num fenômeno de transição, os pa- 
râmetros de que depende o variograma são: 


a— amplitude — a distância a partir do qual as 
amostras se tornam independentes (corresponde à 
noção intuitiva de «zona de influência» de uma amos- 
tra). 


C— patamar — é a variância à priori da V. R. 
estudada e reflecte a «dispersão» própria do fenóme- 
no mineralizado para distâncias superiores à ampli- 
tude. 


C — efeito de pepita— é a ordenada na origem 
do variograma e resulta da sobreposição de microes- 
truturas e de erros de amostragem. 


4.3 — Anisotropias 


Visto que o variograma depende da direcção em 
que é calculado, podem surgir andamentos diferentes 
conforme a direcção. E este aspecto do variograma 
como revelador de anisotropia que dá a informação 
sobre a variação das estruturas presentes com a 
direcção. 

Na Fig. 6 está esquematizada uma anisotropia 
que se verifica numa jazida de fosfatos. 

Num jazigo estratiforme é de esperar que o va- 
riograma calculado no sentido da estratificação seja 
muito diferente (maior amplitude) do que o vario- 
grama calculado na direcção perpendicular em que 
a independência entre amostras estará ligada à es- 
pessura média dos estratos. 

No caso dos parametros C e C (bem como o 
andamento do ramo ascendente) se manterem e só 
variar a amplitude a, trata-se de uma anisotropia 
geométrica que é possível «desfazer» através de uma 
mudança adequada da escala de distâncias que «re- 
produza» o esquema isótropo. 
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Fig. 6 


4.4 — Comportamento na origem 


Para além do caso em que o variograma se reduz 
a flutuações em torno de um patamar (o que signi- 
fica que, à escala de trabalho, não foi possível en- 
contrar a estrutura de relações entre a V. R.), há 
a considerar dois tipos fundamentais de comporta- 
mento do variograma na vizinhança da origem, que 
estão ilustrados na figura seguinte: 


v(h) 


C 


No caso A, diz-se que há comportamento para- 
bólico na origem e verifica-se que a covariância, me- 
dida pela grandeza K(h,) para uma certa distância 
h, é superior à que se observa, para a mesma dis- 
tância, no caso B (comportamento linear). 

Então, um variograma com comportamento pa- 
rabólico corresponde a uma variável mais regular 


no espaço do que um variograma com comporta- 
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mento linear. A covariância dissipa-se mais rapida- 
mente no caso B para as pequenas distâncias, embora 
a amplitude possa ser a mesma. O caso A corres- 
ponde a uma variável como as espessuras de uma 
camada sedimentar, enquanto B corresponde aos teo- 
res de um elemento num suporte constante. 


4.5 — Construção do variograma experimental 


Perante os dados experimentais de que se dis- 
põe, a primeira tarefa que é necessário sempre efec- 
tuar é uma crítica cuidadosa do modo como a amos- 
tragem foi efectuada. Esta crítica exige uma defini- 
ção correcta da V. R. que se pretende estudar, o seu 
suporte e a sua distribuição «macroscópica» no es- 
paço (possibilidade de «viés» por concentração de 
dados numa certa zona, exame de valores aberrantes, 
homogeneidade dos dados (*)). 

Um variograma só tem sentido quando calculado 
sobre variáveis de suporte constante. As variáveis 
devem ser aditivas para que valores médios se pos- 
sam calcular (teores em espessuras variáveis são 
transformados em acumulações, por exemplo). 

No caso das amostras estarem alinhadas e se- 
gundo uma malha regular, o variograma calcula-se 
simplesmente pela fórmula 


1 


H(h)= N' = 


Yes, += 6) | 


onde N' é o número de pares de dados distanciados 
de h, 

Pode suceder que um determinado nó da malha 
não tenha sido amostrado e então o variograma não 


considera os pares de pontos que «liguem» com esse * 


ponto (Fig. 7). 


DADO INEXISTENTE 


3 10 


2 “ 5 6 7 8 E 


Fig. 1 


(1) —> 7 pares 
(2) — 6 pares 
(3) —> 5 pares 
(4) —> 5 pares 
(5) —» 4 pares 
(6) —> 3 pares 
(7) —» 3 pares 
(8) —> 2 pares 
(BD) —> 1 par 


JU] 


Os variogramas calculados em cada linha são 
agrupados em linhas paralelas por média ponderada 
pelo número de pares de pontos que entram em cada 
linha. 


No caso dos dados estarem alinhados mas a ma- 
lha de amostragem não ser regular, definem-se clas- 
ses de distância e(h) e tomar-se-ão todos os pares 
de pontos distantes de h *+e(h) para a estimação 
ds h. e(h) será menor para as pequenas distâncias. 

No caso de dados“aleatoriamente distribuídos no 
plano, faz-se um agrupamento por classes de ângu- 
los. Para a construção do variograma na direcção 6 
e para o passo h, cada valor z(h) é associado com 
todos os que caiem, para a distância h => e(h), no 
ângulo O + 8(0). 


Evidentemente que a prática ditará as classes 
de distância e de ângulos a considerar em cada caso 
concreto. 

No caso extremo em que os processos preceden- 
tes não se adaptam a uma configuração extrema- 
mente irregular de dados, pode recorrer-se à parti- 
ção do espaço em blocos regulares, agrupando toda 
a informação que cai dentro do bloco e calculado o 
variograma sobre as médias ponderadas de todas 
as amostras interiores a cada bloco. 

Todos estes processos de agrupamento de dados 
«regularizam» o variograma experimental e é por 
vezes necessário reconstruir o variograma real. 

O variograma experimental obtido é sempre uma 
curva irregular com flutuações que crescem com h 
(menor significado estatístico, visto que as médias 
são calculadas com menos pares de pontos). Na prá- 
tica, só se tomará um quarto da dimensão total do 
campo, para o ajustamento do modelo teórico. 


4.6 — Ajustamento de um modelo teórico 


Os modelos teóricos de que se dispõe na maioria 
dos casos práticos reduzem-se aos seguintes tipos: 


a) esquemas com patamar 
— esquema esférico 


; dh hs 
C EO h<a 


» h>a 


— esquema exponencial 


| = |h>o 
“h='C|L=ê 1 E 


— esquema de Gauss 


70) = 


alii 
“fua h>o 


a | 


(*) As amostras devem provir de zonas homogéneas da jazida. Por exemplo, não se tratará simultaneamente a parte 


oxidada e os sulfuretos de um jazigo de cobre. 
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b) esquemas sem patamar 
— esquema linear 


yY(h) =p h 


dá 


esquema logarítmico (de Wijs) 


y(h) = A log h h5a 


O variograma mais utilizado na prática é o es- 
quema esférico e o seu ajustamento faz-se em geral 
por um processo gráfico (desaconselham-se os pro- 
cessos automáticos do tipo mínimos quadrados). De- 
pois de obtido o efeito de pepita por extrapolação 
para o eixo das ordenadas, tira-se uma tangente na 
origem à curva experimental. Esta intercepta o pa- 
tamar (que deve ser comparado com a variância à 
priori e ajustado na prática) numa abcissa a' que é 
24 da amplitude a. 


O esquema fica assim definido conhecendo C€, 
C ea. 
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Introdução à estimação geoestatística 


LEOPOLDO PARREIRA CORTEZ (*) 
RESUMO 


O problema da estimação geovestatística consiste na inferência de estimadores de certas características 
globais ou médias de um, determinado domínio (por exemplo, a reserva total ou o teor médio de um jazigo 
mineral), com base na informação obtida através de uma amostragem parcelar, e na quantificação da pre- 
cisão dessa estimação. Esta quantificação é feita através da variância de estimação, e utiliza o variograma 
das variáveis regionalizadas que definem as caracteristicas a estimar, como instrumento de trabalho. 

Neste capítulo fornecem-se as funções teóricas do cálculo das variâncias de estimação e indicam-se 
alguns processos práticos para a realização desse cálculo. A finalizar, apresentam-se duas aplicações con- 
cretas da estimação geoestatística, respectivamente à optimização do reconhecimento de um jazigo filoneano 
sub-vertical por galerias de nível amostradas por rouços, e ao cálculo global da reserva em minério e em 
metal e do teor médio de um jazigo estratiforme amostradn por uma malha regular de sondagens. 


ABSTRACT 


Geostatistical estimation deals with inference of estimators for some global or mean characteristics of 
a certain space based on discontinuous sampling, and with the precision of this estimation, which is quan- 
tified by the estimation variance; the fundamental function in error computations is the variogram. 

Theoretical basis for calculation of estimation variances are referred, and some practical processes 
are proposed. In the last part, two examples of applications are given, respectively to the optimization of 
the sampling of a vein orebody recognized by horizontal galleries, and to the global ore reserve and mean 
grade estimations of a stratiform orebody recognized by a regular pattern of drill holes. 


PLANO: 
1 — Posição do problema. 
11 — O problema da estimação em Geoestatística 


1.2 — Cálculo do estimador 
1.3 — Cálculo do erro de estimação 


Z —Variância de estimação: definição. 


2.1 — Probabilizacão do erro de estimação 

2.2 — Elaboração do estimador 

2.3 — Estimação de uma média por outra 

24 — O variograma como «ferramenta» de estimação». 


3 — Variância de dispersão, 
3.1 — Definição 
3.2 — Relação de Krige 


3.3 — Covariância de duas amostras 
3.4 — Cálculo da dispersão experimental 


(*) Investigador Assistente do I. N. 1, C no C V. R. M. U. L. Eq. a Assistente do Instituto Superior Técnic. 
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4 — Prática da estimação geoestatística. 


4.1 — Generalidades 

4.2 — Cálculo numérico dos variogramas médios 

4.3 — Funções auxiliares 

4.4 — Princípios de aproximação para a composição de variâncias de estimação elementares 
4.5 — O problema do suporte: regularização e acumulação 


5 — Aplicações. 


5.1 —- Aplicação à optimização de um reconhecimento 
5.2 — Aplicação ao cálculo de reservas 


6 — Bibliografia. 


1 — POSIÇÃO DO PROBLEMA 


1.1 — O problema da estimação em Geoestatística 


O problema da estimação em Geoestatística deriva da necessidade de inferior estimadores para os 
valores de variáveis regionalizadas (V. R.) em determinados domínios (os seus campos de aplicação), a partir 
de dados experimentais medidos sobre suportes que não abrangem totalmente tais domínios, ou que podem 
inclusivamente situar-se exteriormente a eles. Trata-se, pois, do problema clássico da inferência estatística, 
que aparece sempre que se pretende estudar uma população através de uma amostragem. 

Como é evidente e conhecido, estas estimações são sempre afectadas de um erro, que é necessário cal- 
cular para que se possa avaliar a precisão da estimação. Este erro é, normalmente, calculado através da 
variância, e pode ser fornecido sob a forma, por exemplo, de um intervalo de confiança válido para um certo 
nível de probabilidade. 

Note-se, finalmente, que a inferência estatística é possível, como se referiu já, devido às hipóteses de 
estacionaridade assumidas para as V. R. 


1.2 — Cálculo do estimador 


Como se sabe, a estatística clássica demonstra que o melhor estimador de uma população de variá- 
veis independentes reconhecida através de um conjunto de amostras é a média (ponderada pelas respectivas 
zonas de influência) dos valores amostrados. 

Mas, como já se viu em capítulo anterior, a estatística clássica não toma em consideração a estrutura 
subjacente ao fenómeno, traduzida pelas correlações entre as variáveis, mas somente os valores amostrados 
independentemente das suas posições relativas, pelo que às estimações realizadas com base na regionali- 
zação da variável correspondem sempre erros menores do que aqueles que a estatística clássica apresenta. 
Exceptua-se, evidentemente, o caso de distribuições puramente aleatórias (efeito de pepita puro), em que os 
resultados obtidos são idênticos por não existir correlação entre as amostras. 

Assim, sucede naturalmente que o melhor estimador de um conjunto de valores de uma V. R. não 
seja uma simples média, geometricamente ponderada, dos valores amostrados na zona a estimar, mas uma 
combinação desses valores que tome em conta as respectivas posições relativas e, portanto, a ligação estru- 
tural existente entre as amostras, podendo eventualmente tomar também em consideração valores de amos- 
tras exteriores âquela zona, desde que ainda estejam correlacionados com ela, 

Este problema de determinar o melhor estimador linear possível do valor de uma V. R. num certo 
campo, com base em toda a informação disponível, foi resolvido por G. MATHERON [(13) ou (14)], que o 
designou por Krigagem, e será tratado noutro capítulo. 


1.3 — Cálculo do erro de estimação 


O erro de estimação, como já se referiu, é normalmente calculado através da variância da distribui- 
ção dos valores amostrados, que designaremos por variância «à priori». 
Esta, como se sabe, define-se pela expressão: 


= E|Z z | 1 É = 
tesfu- = A (n-0) 
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em que Z, são os N valores amostrados, e Z o seu valor médio dado por: 
7 = — 
“0 N 


no caso de a função de ponderação ser constante. 
Também é conhecido que a variância do estimador do valor médio da população é dado por: 


A partir destes valores podem-se calcular, por exemplo, os intervalos de confiança para vários níveis 
de probabilidade. Assim, sabe-se que o verdadeiro valor médio da população tem uma probabilidade de 95 % 


de estar contido no intervalo [Z*—2 6. Z* +29 07. | , que é o nível geralmente utilizado como fornecendo 
uma precisão razoável. 

método da Estatística clássica atrás desenvolvido pressupõe, como se mencionou, que as amostras 
são independentes entre si, o que não é o caso das V. R. Para estas, o método de cálculo dos erros de estima- 
ção tem de ser outro, baseado na teoria das funções aleatórias, que se expõe a seguir. 

Para tal, convém fazer uma primeira distinção entre dois tipos de varincias: a variância de esti- 
mação (notação 57), que caracteriza a qualidade de uma estimação, e a variância de dispersão (notação D:), 
que caracteriza a variabilidade espacial de uma quantidade (de um teor, por exemplo). 


2 — VARIÂNCIA DE ESTIMAÇÃO: DEFINIÇÃO 
2.1 — Probabilização do erro de estimação 


Seja z o verdadeiro valor desconhecido de uma característica de uma determinada zona do campo 
onde ela é definida, que se pretende calcular pelo estimador z* (de valor conhecido através de uma Kriga- 
gem, ou pela média de valores amostrados, ou mesmo por uma única amostra, por exemplo). O erro que se 
comete na estimação de z por z* é (z-z*), e poderá ser caracterizado, em termos probabilísticos, pela sua 
variância. Vejamos como: 

O método probabilístico usado na Geostatística interpreta z, z* e (z— z*) como realizações parti- 
culares de funções aleatórias (F. A.) Z,Z* e (Z-Z*), isto é, como V. R. 

Esta interpretação é possível graças às hipóteses de estacionaridade já referidas, nomeadamente à 
hipótese intrínseca que permite admitir que as características estatísticas de um dado domínio permanecem 
as mesmas, em média, quando este se desloca por translação através do espaço. 

Assim, se deslocarmos no espaço a zona (um painel mineiro, por exemplo) e as amostras que servi- 
ram para a construção do estimador z*, ligados um ao outro, a cada implantação do painel corresponde 
um valor para o seu teor verdadeiro z, e um só valor do estimador z*. 

Quer dizer, o erro (Z-Z*) é uma variável aleatória que pode ser caracterizada pelos seus dois pri- 
meiros momentos: 


— a média: m = E[Z-Z*] 


— a variância de estimação: o, = E ( [ Z— Z* ) — mº 


Esta variância de estimação de Z por Z* é portanto, a variância do erro (Z-Z*), em que cada valor 
numérico é tomado numa implantação particular do conjunto ligado painel-amostras. 

Devido à hipótese intrínseca, esta variância permanece constante se for medida numa outra zona do 
mesmo domínio, mas varia se for calculada para outra posição das amostras na zona (painel). 

Note-se que sendo m -£O, a estimação é enviezada. Sem = O, a variância de estimação será dada por: 
o =E [Z2-2Z*]2. 


2.2 — Elaboração do estimador 
O estimador Z* é sempre, de um modo ou de outro, uma função de informação I disponível através 


da amostragem realizada. Suponha-se, por exemplo, que I seja um conjunto discreto de N amostras z(x,), 
i=1aN. | 
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Então, será: Z* = f[Z(x,), Z(x,), ...Z(X,), ...Z(Xy)] 
Esta função de IN variáveis não pode ser qualquer. Tem, nomeadamente, de respeitar as duas condi- 
ções seguintes: 


a) deve verificar a condição de não enviezamento: E[Z-Z*| = O 


b) deve ser tal que se possa calcular a variância de estimação, quer dizer, os termos do desenvol- 
vimento: 


E[Z-2*|*= B[Z”] + E[2**]|-2E[Z2*] 
É devido à segunda condição que se utilizam geralmente estimadores lineares, do tipo: 


O = à Z(x) 


Um estimador de certo tipo diz-se óptimo se minimiza a variância de estimação, por exemplo, o 
estimador de Krigagem é um estimador linear, não enviezado e óptimo, 


2.3 — Estimação de uma média por outra 


Para fixar ideias, estudemos o caso particular da estimação do teor médio z, num domínio v através 
do teor médio z, calculado sobre outro domínio v', estando ambos os domínios compreendidos num campo 
homogéneo. 


Então, admitindo uma distribuição uniforme da V. R., será por definição: 


Z, = | Z (x) dx 


Y 
. os 


* 


| 
v = | ” Z (x) dx 


Onde z(x) designa a V. R. pontual, e seja Z(x) a F, A. correspondente, suposta estacionária de 2.º 
ordem e com variograma y(h). 


Note-se que os domínios v e v' podem ser definidos no espaço a 1, 2 ou 3 dimensões; por exemplo, no 
espaço a 3 dimensões, dx será dx, dx, dx,. 


Admitindo a hipótese estacionária, não há enviezamento, pois E[Z | = E[Z,,|). Portanto, a variância 


do erro (Z,-Z,) será: o = E[Z,-Z,]º = E[Z5] + E[Z,º] — 2E[Z,2,). Mas a variância de Z(v) é dada 


9 1 * ” 1 - - | 
E | Z | = 4 | | É [7 (x) . Z (y) | dx dy = çã | Í C (x-y) dx dy 


por: 


. E." v 


em que: C(x—-y) — E[Z(x).Z(y)] é a covariância da V.A.Z quando x e y descrevem, independentemente 
um do outro, o domínio v. 


De igual modo: 


ER ca q 1 
| 2] = Es | . | e[Z0.26) | dx dy E) 


percorrendo agora x e y, independentemente um do outro, o domínio v'. 
A covariância de Z com Z,, será: 


| 1 * - Ed í - - 
E| 2, 2] — =] | | 09:26) | dx dy = ep | | C (x-y) dx dy 


“ v. v 


| | C (x-y) dx dy 


A v'd v 
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onde x e y percorrem, um, o domínio v, e o outro, o domínio v'. 


Logo, substituindo valores: 


* 


' Tr Ww 1 | 2 
E Z, — 2] A J C (x-y) dx dy + + |, [ C (x-y) dx dy — =| | C (x-y) dx dy 
E Lori v Po o W' va v' 


Mas, conforme já se referiu, o variograma pode-se escrever: 


y (h) = C(O) — C(h) 


Donde, fazendo |x = il =h, e substituindo C(h) por C(O) = y(h), como C(O) desaparece, vem a 
expressão final: 


4 É ) a 1 cs 1 cm | 
GE) L, — Z,| = J [ 7 (x-y) dx p=] )J 7 (x-y) dx gy), | 7 (x-y) dx dy 
Ver va v' v JdoYw Y v w' v' 


Ou, em notação simbólica: 


Pois) y(s 7 (vv) 


O símbolo y(v,v'), por exemplo, designa o valor médio do variograma y(x—-y) quando os dois pontos 
de apoio x e y descrevem indepedentemente, um o domínio v e o outro o domínio v'. 

Esta expressão tem a vantagem de permanecer válida para qualquer F. A, intrínseca, mesmo se a 
covariância C(x—y) não existir. Todavia, se se preferir trabalhar com covariâncias (e se elas existirem) 


pode-se usar a notação simbólica: 


AC, rT(w,v)—-2T(v,v) 


Estas expressões podem estender-se a domínios v e v' de qualquer geometria e não forçosamente 
compactos ou contínuos: por exemplo, o domínio v a estimar pode ser constituído por dois painéis distintos, 
v=v, +v,; os domínios v e v' podem ser disjuntos (total ou parcialmente) ou estar um incluído no outro; 


o domínio v' pode ser constituído por várias amostras discretas, etc. 
A leitura das expressões de e mostra 4 factores essenciais e intuitivos de que depende a qualidade 
E 


de qualquer estimação de v por v': 


a) a geometria do domínio a estimar: termo /(v,v) 


b) a geometria do domínio estimante: termo y(v',v') 
c) o grau de regularidade do fenómeno estudado: utilização da característica estrutural variograma 


d) a geometria relativa (distâncias) de v e v': termo Y(v, v'). 


O modo de cálculo dos valores médios do variograma será tratado à frente. No caso, atrás referido, 
de v ser constituído por N amostras discretas x,, e z,, ser a média aritmética desses valores, isto é: 


1 N 
Dé 3 z (x,) 
1 


N 
Será: Tt) = — BD V(x— x, ) dx 
N.viã dy 


quis l 
Tv) == f fi rixn azar 
N 


Í N 
2 a | (x — Xj) 


Es o JA 


“e] 
a, 
ci 
a 
e 
| 
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A fórmula atrás deduzida também se aplica, como é evidente, à variância de estimação do teor 
médio z de um painel v por uma combinação linear z das informações disponíveis. 


Por exemplo, se se dispõe de N informações S, de teores médios z(x,), e sendo à, o ponderador asso- 
ciado à informação S., o estimador é, como já se viu: 


Esta expressão é geral quaisquer que sejam as geometrias do painel v e das informações S, e quais- 
quer que sejam os ponderadores à,. Contudo, o não enviezamento deve ser assegurado: E[Z — Z*| = 0, 


para o que basta impor a condição: E à, = 1, como se verá no capítulo seguinte sobre Krigagem. 
] 


24 — O variograma como «ferramenta» de estimação 


As expressões atrás deduzidas têm a vantagem fundamental de necessitarem apenas do cálculo de 
variogramas médios para se conhecer o erro de qualquer estimação. Isto é, realizada uma análise estrutural 
do campo de aplicação da V. R. em estudo, e definido o esquema teórico (o modelo) que melhor interpreta 
o variograma experimental da regionalização (observando todas as condições já mencionadas no capítulo 
respectivo), podem-se avaliar os erros inerentes a qualquer estimação feita dentro do domínio de validade 
desse modelo. 

Esta propriedade não é de estranhar visto que, como se sabe, o variograma (pela sua própria defi- 
nição: 2 y(h)=E[Z(x+h)-Z(x) |2) pode ser interpretado como a variância do erro que se comete ao esti- 
mar um teor desconhecido (em (x+h)) por meio de um dado z(x) (em x). Como em estatística se carac- 
teriza a amplitude de um erro aleatório pela sua variância, pode-se considerar 2 y(h) como a variância de 
estimação de z(x+h) por z(x) (ou vice-versa). 


3 — VARIÂNCIA DE DISPERSÃO 


3.1 — Definição 


I4- 


Considere-se um conjunto discreto de N valores Zz(x,), com i=1 a N. Seja: Z- 


z (x,) 0 valor 
médio deste conjunto. 


Z|m 
n 

| 
—. 


Os N valores z(x.) dispersam-se em torno da média z, podendo esta dispersão ser caracterizada pela 


expressão: 
N 2 
> | Zz(x,) — z | 
= ; 


Tendo o conjunto dos N valores sido tornado aleatório, pode-se definir a variável aleatória d: cuja 
esperança é, por definição, a variância de dispersão do conjunto discreto das N V. A. 4(x,), i=1 a N: 


3 as 


dd Ra 


1 á 

" : 72 

Dº - a E| E Z(x)— Z] | 
i 

A dispersão experimental dz representa então um estimador desta variância de dispersão De. 


Esta definição pode ser generalizada a conjuntos contínuos do seguinte modo: seja um painel qual- 
quer de volume V centrado no ponto x. O teor médio deste painel será: 


- 


| 
Ze (X) = V | z (x + u) du 
J wv 
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Considere-se, no interior deste painel V, um bloco (unidade de produção) de volume v centrado no 
ponto w, e de teor médio: 


1 a 
Zy (3) = o) z (y + w) dw 
W 


A cada posição y da unidade v no interior do painel V, corresponderá um valor particular z (y); a 
dispersão destes valores em torno da sua média z (x) pode ser caracterizada por: 


sa 


1 
()= | My)? dy 
dv 


O conjunto dos valores z (y) e z (x) tendo sido tornado aleatório, então sº(x) aparece como uma 
realização particular da V. A. S:(x), dada por: 


Sº (x =). [Z,(y)—Zy (0) 7 dy 


A esperança (valor médio) desta V. A. S:(x) é, por definição, a variância de dispersão do volume v 
no volume V: 


* 


1 
D* (v/V) = E | v | [Z2,(y) — Zy (x) PP dy | 
Jv | 


A hipótese da estacionaridade da F. A, Z(x) implica que E[S:(x) | não depende da posição particular 
x do painel V, e portanto a expressão de D:(v/V) é geral e pode ser calculada do seguinte modo: consi- 
dera-se o volume V como a reunião de N volumes v, disjuntos, todos iguais entre si e de valor v (ou seja, 
cada um dos V, deduz-se de v por translação). Seja Z(x,) o teor médio de Z(x) em VW e Z(V) o teor médio 
de volume V. Então, de acordo com o que se referiu, vem: 


l 


1 
DG =D= q) »[Z(x)—Z (V)P N ds É — Z (V)P 
] Í 


Ou seja, é o valor médio das variâncias de estimação de v, a V. Donde, notando que V = LE v, = Nv, 
o desenvolvimento da expressão vem: ! 


D* (v/V) "2! E ([Z(O)P)HE(IZ(WP)-2E[Z(x).Z(V]]= 
] 


1 1 Fr 

No ESA E» C (x,- nn gi. CG (x-y) dx dy — 
2 

o Ng El Vi NE fis Pas, g= |, f da 
LN 


C (x-y) dx dy 


=J,), 
Ou seja, em notação simbólica: 
DV) == E (ro) CT (VV) 


Ou ainda, utilizando o variograma de F. A. pontual Z(x) em vez da covariância: 


MO À IN 7 (9) day — fo fores )dxdy=7(V,W)— 7 (v,v) 
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Isto é, a expressão da variância de dispersão depende apenas do variograma e das geometrias respec- 
tivas de ve V. 


3.2 — Relação de Krige 
Suponha-se o caso da variância de dispersão de v em V, quando o volume v está reduzido a um ponto. 
Então, a variância D:(0/V) da variável pontual Z(x) é dada por: 
D*(0/V) = E[Z(x) —- Z(V)]” 
Ou seja, desenvolvendo: 


D(OVW)=E([Z()P)J+E ([Z(VW)P)—-2E[Z(M).Z(V)] = 


- . ç 


5) - 
= | | | | C (x-v) dx dy — V | C (x-v) dy 
Val o Y 


“ 


Atendendo a que: 


e | C (x-v) dx dy 
Vdvdvy 


Vem finalmente (passando para variogramas): 


DON) == | t(xy)dray=7 (VV) 
Vo Vv 


Isto é, o valor médio do variograma quando os dois pontos de apoio descrevem, independentemente 
um do outro, o volume V, 


E evidente que será também: 
Dº (O/V)=7 (vv) 
Donde se pode escrever a expressão da variância de dispersão de v em V: 
Dº (v/V) =D? (0/V) — Dº (O/v) 


Designe-se por G o volume de um jazigo composto por vários painéis de volume V, dividido, por 
exemplo, em unidades de volume v, de tal modo que vcVcG. Neste caso, será: 


D> (v/G6) = Dº (0/6) — Dº (O/v) 
Subtraindo a expressão anterior desta: 
Dº (v/G) — Dº (v/V) =Dº (0/6) — Dº(O/V) =Dº (V/G) 
Ou seja, a relação de aditividade: 
Dº (v/G) = Dº (v/V) + Dº (V/G) 


Isto é, a dispersão da unidade v no jazigo G é igual à soma da dispersão de v no painel V com a 
dispersão de V em G. 

Esta relação é conhecida como relação de Krige, por ter sido experimentalmente verificada por este 
técnico nas minas de ouro sul-africanas, bastante tempo antes de ter sido estabelecida teoricamente. 

Verifica-se assim que, qualquer que seja a posição do painel V no jazigo G, as suas unidades internas 
possuem a mesma variância de dispersão (resulta de D'(v/V) =D? (v/G) —- D'(V/G)). Isto é, esta não 
depende da localização particular do painel, mas única e exclusivamente do seu tamanho e forma. 

Diz-se então que o espaço mineralizado possui uma lei de dispersão intrínseca, ou seja, que o grau de 
descontinuidade da mineralização é sempre o mesmo (V. hipótese intrínseca). 
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O 


3.3 — Covariância de duas amostras 


Sejam v e v' duas amostras de distância e disposição mútua fixas no campo YV. 


Define-se, de forma análoga às anteriores, a sua covariância como: 


“ 


º l | is E 
D (v,v/V)= E] | Í 7 (x-y) dx dy — RR [ [ | 7 (xy) dx dy =7(V,W — 7 (v,v) 


Assim, a variância de estimação pode ser expressa em função das variâncias e covariâncias de v e v' 
em V por: 


6,=Dº (v/V) + Dº(v/V)— 2D (v, v'/V) 
3.4 — Cálculo da dispersão experimental 


Se se dispõe de N dados experimentais z (x,), i=1 a N, repartidos uniformemente no campo V, a 
dispersão experimental: 


é um estimador do valor teórico da dispersão Dº(v/V). 


Convém notar finalmente que os resultados teoricamente obtidos traduzem dois fenómenos experimen- 
tais bem conhecidos dos mineiros: 


a) A dispersão de um conjunto de dados amostrados num volume V em torno do seu valor médio, 
cresce com a dimensão de V; tal é a consequência lógica da existência de correlação espacial 
entre as amostras: quanto mais próximas estão, mais os seus valores são, em média, próximos 
uns dos outros. 

b) A dispersão de um conjunto de dados de volume v, amostrados num volume V fixo, em torno 
do seu valor médio, diminui quando v cresce; é o fenómeno da regularização: a variabilidade dos 
teores médios dos blocos de exploração é atenuada em relação à dos teores de «carottes» de sonda- 
gens, por exemplo, 


4 — PRÁTICA DA ESTIMAÇÃO GEOESTATÍSTICA 
4.1 — Generalidades 


As expressões atrás deduzidas são, normalmente, de difícil aplicação prática. Note-se que, por exemplo, 
o cálculo do variograma médio: 


= 1 o + 
7 (v,v)=—— J | rtp dxay 


implica o cálculo de um integral sextuplo (num espaço a 3 dimensões), cujo cálculo analítico directo é 
complexo. 


Duas soluções práticas se apresentam: 


a) Cálculo numérico (em ordenador) dos valores médios dos variogramas. 


b) Decompor a resolução analítica dos integrais sextuplos em etapas sucessivas, algumas das quais 
são previamente resolvidas e tabeladas uma vez por todas. Estas etapas intermédias correspondem 
à definição de funções auxiliares e à utilização de ábacos, e a sua interligação implica a admissão 
de princípios de aproximação para a composição dos vários termos. 
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Note-se ainda que o problema do suporte das amostras, como se referiu em capítulo anterior, tem 
implicações no cálculo dos variogramas. Com efeito, as expressões deduzidas são válidas para variogramas 
pontuais (com suporte de volume nulo ou que possa ser tomado como desprezável em face do volume a 
estimar), enquanto que, na prática, não se tem nunca acesso a valores pontuais z(x,), mas sim a dados 
regularizados sobre um certo suporte de dimensões, por vezes, consideráveis. 

É, então, necessário deduzir os variogramas pontuais a partir dos variogramas regularizados. 

Este problema da regularização será tratado adiante. 


4.2 — Cálculo numérico dos variogramas médios 


É a solução normalmente mais rápida se se dispõe de um ordenador, ou mesmo de uma pequena 
calculadora programável. 

Implanta-se uma malha regular (ou aleatória, ou aleatória estratificada) x,, com i=1 a N, no volume 
v, uma outra (idem) Y, com j=1l a N', no volume v, e assimila-se o integral sextuplo a uma dupla soma 
discreta: 


N N' 


= - 1 
Vir) e Pim ) = N N' bd » 7 (x) 
Ev =| j= 


Conhecidas as coordenadas x, e y,€o modelo y(h), o ordenador calcula automaticamente este esti- 
mador do valor médio do variograma. 

Um exemplo de programa deste tipo é dado em M. DAVID (5), existindo também outros na biblio- 
teca de programas do C. V. R. M. U. L. 

O erro deste estimador estã dependente, como é óbvio, da densidade de discretização no interior de v 
e v, diminuindo quando N e Nº aumentam. É essencial que este erro não seja tal que mascare a variabili- 


dade estudada. Assim, para uma variância de dispersão D* V/V)=Y(V, W)-—y(v, v) os estimadores numé- 
ricos y*(V, V) e y*(v, v) devem verificar a relação de ordem teórica Y(V,V) >Y(v,v). Também para uma 


variância de estimação gq? =2y(v, v)—y(v, v)>y(v', v'), os três erros numéricos arriscam, ao acumular-se, a 
[5 


mascarar a variância gq” procurada. 
E 
Na prática, as três regras seguintes devem ser sempre observadas: 
a) Os diferentes estimadores numéricos devem conservar entre si todas as relações de ordem teóricas. 


b) A discretização deve ser a mesma para todas as estimações numéricas y* dos valores y de uma 


mesma fórmula. Assim, por exemplo, os três valores y da expressão q =2y(v,v)-y(v, V)>y(w', v') 
E 
devem ser estimados com a mesma discretização (N pontos em v e N' pontos em v'). Esta regra 
garante a coerência da fórmula: com efeito, a variância é é então calculada como a variância de 
E 
estimação de N pontos de v por N' pontos de v'. 
c) A densidade de discretização depende dos objectivos visados, e será em geral optimizada com a 


experiência. De início, pode-se proceder por iterações, até que o suplemento de discretização já 
não melhore significativamente a precisão da estimação. 


Como regra empírica, e para fins de Krigagem, bastam, em geral, as seguintes discretizações: 
— N=10 pontos, para v a 1 dimensão 
— N=6X6 pontos, para v a 2 dimensões 


—N=4xX4x4 pontos, para v a 3 dimensões. 


No capítulo sobre krigagem, onde o cálculo dos valores médios do variograma é essencial, este ponto 
será retomado e aprofundado, 


4.3 — Funções auxiliares 


Uma função auxiliar é um valor médio y(v,v') correspondente a geometrias v e v' particularmente 
simples mas de uso frequente nos estudos habituais. Normalmente, usam-se 4 funções auxiliares, simboli- 
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“4 Guindastes eléctricos para o cais comercial do Porto de Aveiro: 


— 1 Guindaste de 6 — 121/24 — 20 m 
— 3 Guindastes de 3 — 6t/24 — 20 m” 


PONTES ROLANTES, GUINDASTES E 

APAR. DE ELEVAÇÃO ESPECIAIS Projecto e fabrico 
TURBINAS HIDRÁULICAS —— Fabrico segundo licença de A.C. M. de Vevey, S. 4. 
- TURBINAS A VAPOR ——— Fabrico segundo licença de Brown Boveri, Cie. 
CALDEIRAS À VAPOR ————————— Projecto e fabrico segundo licença de Foster 
Wheeler, Co. 


EQUIPAMENTOS E INSTALAÇÕES 
INDUSTRIAIS 


CONSTRUÇÕES METALOMECANICAS MAGQUE ;..-. 


ALVERCA DO RIBATEJO - PORTUGAL: 


TÉCNICA NV 


PLA 


y 


o que é 


PLA 


Você sabe 
know-how 


TÉCNICA VI 


& f1tj 


Slas Copco 


zadas y, F, qa e H, definidas no espaço a 1 dimensão (as duas primeiras) ou a 2 dimensões (as quatro), 
consoante a dimensão da figura de integração. 


a) Espaço a 1 dimensão: 


A B x = MM 


-—o 0 el= == e 2 caga 


pano [es 


Seja o segmento AB de comprimento L, e y(x) o variograma 


pontual neste espaço. 


— Valor médio de y(x) quando M está fixo em A e M' descreve o segmento AB de comprimento L: 


E O si 


— Valor médio de y(x) quando os 2 pontos M e M' descrevem, independentemente um do outro, o 
segmento AB de comprimento L: 


a eg do ri 1 “IJ, ml, 
| ()=7 (AB,AB) = 5 | | 7 (|x-x' |) dx dx 
“o! 0 o 


Demonstra-se que este integral duplo se reduz aos integrais simples seguintes: 


j 2 E a | 
E je | Se duim di | (L-x).7 (x) dx 
= ao 0 oo 


bj Espaço a 2 dimensões: 


po 


h=— MM 


Seja o rectângulo ABCD de lados (L, 1) no espaço a 2 dimen- 
sões, e seja y(h) o variograma pontual isótropo (ou tomado como 
tal, depois de desfeitas eventuais anisotropias geométricas) neste 
espaço, 


C — Valor médio de y(h) quando M descreve um lado de com- 
D primento 1 e M' descreve o lado oposto 1, distante de L;: 


| E quê A Estes 
a(L,)=7 (RC, BD) =) (=x) 7 (WIZ =) dx 
LH) 


De modo análogo, define-se também a(l,L)=y( AB, CD). Note-se que esta função q não é simétrica: 
a(L, 1) + a(L L). 


— Valor médio de y(h) quando M descreve um lado de comprimento | e M' descreve todo o rectân- 
gulo (L,1): 


as | , O ia 
X(L,)=7 (AC,ABCD) = (BD, ABCD) => [ a (x,1) dx 
Ná 


Esta função também não é simétrica: y(L,1) & y(1, L). 
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— Valor médio de y(h) quando M e M' descrevem, independentemente um do outro, todo o rectân- 
gulo (L,1): 


Rai sc - dia A iii 
F(L,D=F(1,L) = 7 (ABCD, ABCD) =77 | x.Z (x,1) d= f, (L-x). 2 (x,1) dx 
a 0 


— Valor médio de y(h) quando M descreve um lado de comprimento 1 e M' descreve um lado adja- 
cente de comprimento L;: 


H(L,)=H(,L)=7 (AC,AB)=7 (AC,CD) = 7(BD,AB) =Y(BD,CD) = 
1 dq 1 o? 
— Pesa = 2 b = a Mu pos L.lr.F Es 
21 91 (IE. X (Ls 1) 4LiI aldL CL 1) 


Demonstra-se que H(L,1) é também o valor médio de y(h) quando M está fixo num vértice e M 
percorre todo o rectângulo (L, 1): 


H(L,1) = y(A, ABCD) 


c) Equivalências: 


Note-se que basta reduzir a zero um dos lados do rectângulo ABCD para fazer corresponder as fun- 
ções auxiliares a 1 e a 2 dimensões: 


Q(L,0)=7/(A,B) =7(L) x (0,)=7(AC,AC)=F (1) 
L(L,;O)=7(A,AB) =X(L) x (0,)=7(AC,AC)=F (1) 
F(L,0)=7(AB,AB)=F(L) F (0.)=7(AC,AC)=F (|) 


H(L,0)=7 (A,AB) =YZ(L) H (0,)=7(A,AC) =Yx(1 


Observe-se ainda que todas as funções auxiliares a 1 ou 2 dimensões se podem deduzir duma só, 
F(L, 1), por derivação ou anulação de um parâmetro: 


x (l..l ! cd Lº.F (1.,1)) 

(1, ) 2 9 12 (L. o (1,1) 
z l 4 ) 21 1. (1 º (1, )) L(L, ) 21, 91, ! a «1) 
H(L.l) = : a L*,B.F(L,] F(L)=F(L,0 

cocada 4L1 01.9] (Le, dE. E (d.,1)) (L.) = «3 4) 


d) Espaço a 3 dimensões: 


Existem também as correspondentes funções auxiliares q, y, H e F para configurações paralelipipé- 
dicas a 3 dimensões, mas na prática não são utilizadas. 

Empregam-se ábacos para as funções auxiliares definidas sobre paralelipípedo de base quadrada (LX 
X1xX1), dados, por exemplo, em A, JOURNEL (9), em M. DAVID (5), ou em J. SERRA (17), ou então 
usam-se os princípios de aproximação para a composição das variâncias. 


e) Expressões e Abacos das Funções Auxiliares para Esquemas Isótropos 


Em A. JOURNEL (9), M. DAVID (5), J. SERRA (17), F. MUGE et al. (15), entre outros trabalhos 
listados na Bibliografia, apresentam-se expressões e ábacos das funções auxiliares válidas para os esquemas 
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mais vulgares (linear, logaritmico ou de De Wijs, esférico ou de Matheron, parabólico ou de Gauss, exponen- 
cial ou de Formery, etc.), válidos para campos isótropos no espaço a 3 dimensões. 

A título exemplificativo, indicam-se as expressões das funções auxiliares a 1 dimensão para o esquema 
esférico normalizado (efeito de pepita C —O, patamar C=1, alcance a): 


à h 1 h3 

a Yhe[0,a] 
Pos; 2". “a 

1 Yh>a 

3 1 1 148 

4a 8a rs AS 
X (L) = ' 

r 1 

1—-— — Es 

| 5 L Wiki 

1 É. 1 Is 

2 a Ma PR [a 
F(L)= 


ias WI a 


L, 1 
Nos ábacos n.º* 1 a 6 estão tabelados, em função de — e de —, as funções seguintes: 
a a 


Abaco n.º 1: a(L,1) e a(O,1) 
» D.º 2: x(1 1) 


>» n.º 3: H(L,1) 

* n.º 4: F(L,1) 

» n.º 5: F(L,1,1) 

» n.º 6: a(L,1,1) e «(0,1,1) 


Os ábacos n.º 7 a 10 referem-se às variâncias de extensão elementares e são explicados adiante. 
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ABACO N.2 
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ABACO N.3 
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ABACO N. 4 


função F(L,|) 


Po 
CAS 
SPA? 
TAS 
ps) A A? 
|) A Ro 
oRa Ao 
PE 
E vo 
VAGA 
Po 
S9/9 
OG j 
So 
S% 
o" py | 
eh E o 
Sds 
Sto 
EA 
“Fa 
o pP É 
o 8º 
Ss +» 
qo Ea 
SP. 
o” 
Q 
somou 4 io N 
E jo “Doooo [a o o 


5 6 789) 


4 


04 0.5 0.60.70.8091 


0.3 


0.1 


Seg. A. JOURNEL (9) 
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ABACO N. 5 


função F(I,h,h) 


od Do 
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Seg. M. DAVID (5) 
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ABACO N. 6 
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ABACO N. 7 
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Seg. A. JOURNEL (9) 
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ABACO N. 8 
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HT 

+ 

= 

casi 

o] 

- 


AU iai aa a 


oz ii ha K mê 


0.1 0.2. /1 L A ui] | l 5 6 7 8910 


Seg. A. JOURNEL (9) 
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ABACO N. 10 


ho 


0.2 


0.1 
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tj Cólculo das variâncias de estimação elementares. 


A escrita das variâncias de estimação utilizando as funções auxiliares é imediata. Assim, apresentam-se 
alguns exemplos de uso corrente: 


— À 1 dimensão: 


Estimação do segmento L por uma amostra pontual central: 


, A () B 

= - = : — fold EEE na me e e e e 

cp — 2 7(0,AB) — y(AB, AB) — 710) = 2 x má ) — F(L) R——>— L———— 
Estimação do segmento L pela média de duas amostras pontuais extremas: A B 


1 
é = E s(u— PL = — (Lo 
E k > Y 


A 1 B 
je 2 ma o if 
Estimação do segmento L, por um segmento | central: pe il a 
o é 1 
qp=2% — ) = P(L) — F(1) a | | 18 
2 


Estimação do segmento L por dois segmentos | extremos: 


1 | 
e; = 2 y(L) — F(L) — = | o + em] 


— A 2 dimensões: 


Estimação do rectângulo (L,1) por uma amostra pontual central: 


bu 4: 
op =2H [=,5) — F(L,1) 


A B Estimação do rectângulo (L,1) pela média de 4 amostras pontuais nos vértices: 
| [o] | i 
1 f 4 
ce = 2 H(L, 1) — F(L, 1) — — | y(L) + y(1) — T+? | 
Cc L D E EI ria 
| ES Estimação do rectângulo (L,1) pela sua mediana 1: 


SR 
cp =2 Arts — F(L, 1) — F(1) 


À 
nm Estimação do rectângulo (L,1) pela média dos seus dados 1: 


1 1 
dp = & x(L,d) — FI(L,1) — la) pás 5 FO) 
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— À 3 dimensões; 


Estimação do paralelipípedo (L,1,1) pelo eixo L;: 


1 1 
0; = 2H (5) —F'(L,1,1) — F(L) 
2 3! 


L/2 L/2 


Estimação do paralelipípedo (L,1,1) pelo plano mediano (1,1): 


L 
q =2% (ti) - F(L,1,D — F(,1) 


Nos ábacos n.º* 7 a 10 estão tabeladas algumas destas variâncias de estimação elementares. 


44 -— Princípios de aproximação para a composição das variâncias de estimação elementares 


Como já se referiu, admitem-se geralmente um certo número de hipóteses, ou princípios de aproxi- 
mação em Geocestatística, com o objectivo de simplificar os cálculos sem lhes diminuir significativamente 
a precisão, 

Um desses princípios utilizados na estimação geoestatística é o princípio de independência dos erros 
elementares, que permite a sua composição directa. Consiste em reduzir o erro global (Z —Z*) de uma esti- 
mação a uma soma de erros parciais elementares considerados como independentes entre si e cuja variância 
é simples de calcular, 

Assim, seja o erro global: (Z -Z*) = 54 h; (Z,—2Z,*), em que à, é a função de ponderação. 


A variância de estimação global será; 
q =E[Z-2*])'=0 NX E[Z,-2*]) + DOM E[(Z-2*) (2,-2")] 
5 


Se os erros (Z,—-Z* ), (2,-2.*), com i-£j, poderem ser considerados como independentes, todos os ter- 
mos rectangulares (covariâncias) se anulam, e fica somente: 


N N 

” ; = 2 

cg =0 ME[Z-2*]) =L.& o 
i=1 i=1 


As configurações das estimações elementares de Z, por Z * serão escolhidas suficientemente simples para 
que se possam calcular directamente as variâncias EA por meio de fórmulas gerais ou com o auxílio das fun- 
ções auxiliares ou de ábacos. 

Pode-se avaliar o grau de aproximação da hipótese de independência calculando os N(N-1) termos 
rectangulares tais que: 


E[(Z-2º) (2-2) ] =E[2, 2,)]-B(Z, 2º] -B[Z, 2º] +B[Z8 28] (com 14) 


pois cada um dos termos deste desenvolvimento exprime-se em função dos valores médios do variograma 
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Distinguem-se dois tipos de composição de erros elementares: 


— Os que provêm de uma mesma configuração (composição directa) 
— OS que provêm de configurações diferentes (composição de termos). 


a) Composição directa de erros elementares: 


Seja, no espaço a 1 dimensão, uma linha L reconhecida por N pontos regularmente espaçados de 1. 


O 
— Dispositivo centrado: MS O O 
a 


A linha L é constituida por N segmentos 1 reconhecidos cada um por uma informação (amostra) 
central: 


A variância da estimação elementar do teor médio Z, de cada segmento 1 pela informação central 
Z(x,) é dada pela expressão geral: 


: ; f AS 
o; Ed r (tl x )—-y(x, x) —y (1, 1) E a 2 =P (1) (v,) 


O teor médio Z, de L será dado pela média dos teores Z, em cada segmento 1.: 


1 N 
E =— já XxX) 
N i=1 
E o seu estimador por: 
e 1 N 
Z = — Z(x.) 
| N b. Ê 


Ss 
A variância da estimação global de Z, pela média aritmética das N amostras Z(x,) vem: 


x 


o; = E[Z,-Z, |] = 


o MEl-2(]) 


Supondo os erros de estimação elementeares (Z, —Z(x,)), (Z,—Z(x,)) como independentes, y 1/3, 
vem: 


2 mo O E[Z,-2 ms É 
CE = — -Z(x)]'= — 0 
E. Nº | | Hi Í | N “| 


— Dispositivo fechado: 


A mesma linha L pode também ser estimada por (N+1) amostras em dispositivo fechado (segmen- 
tos amostrados nos extremos). 


Neste caso, o teor médio Z, de cada segmento | é reconhecido pela média dos teores dos seus extremos: 
* 1 
4 = = [Z2(x,) + Z(X,) | 


Então, o estimador do teor médio Z, de L, vem: 


“a” as 1 Z(x,) Z(XK 1) 


Re E Sa 
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Fitas auto-adesivas Tesa 
para a indústria eléctrica 
e electrônica. 


Ja nar 
AA” 
Ds ” 
N EN 
A 


Se 


1, 


Ea 
OA, 


RR 


in 
= ia 


As fitas auto-adesivas Tesa técnicas referentes a eléctrica e electrónica. Uma 


para a indústria eléctrica 

e electrónica oferecem todas 
as vantagens das fitas, que 
foram estudadas para um fim 
específico. E nas mais diversas 
finalidades: para isolar, ligar, 
fixar e referenciar. Elas 
consoante a-sua aplicação 
satisfazem as especificações 


isolamento térmico, 
resistência em superfícies de 
verniz, inalterabilidade às 
variações atmosféricas e ao 


ataque de lexívias ou ácidos. 
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TÉCNICA VIII 


Mas agora, dois segmentos contíguos 1 e 1, têm em comum uma mesma amostra, pelo que já não é 
possível considerar os seus erros de estimação elementares como independentes, 

Todavia, demonstra-se que é possível, com boa aproximação, assimilar um dispositivo fechado com 
(N+1) pontos a um dispositivo centrado com N pontos, desde que N > 5. Esta regra prática estende-se a 
dispositivos fechados a 2 ou 3 dimensões, 


” cas i = e] 
— Dispositivo aleatório: ps L 
pu | gem 
o qua 


Considere-se agora a mesma linha L reconhecida por N pontos em dispositivo centrado mas de malha 
irregular (1-1), sendo cada um dos segmentos 1, estimado pela sua amostra central Z(x,). 


| - +. 8 
Então, o teor médio em L será: Z, = E 1 l, Z, 
j=1 
+ 1 N 
E o seu estimador: Z = = Li Z(x,) 


A composição directa dos erros de estimação elementares (Z,-Z(x)) dá: 
É RE 2 E E É.ã 
cp — É[4,—4, | = EFE +. À B[4-4(x,)]' = “Tr lo 


E, contudo, necessário que a malha seja aleatória (pura ou estratificada), isto é, não esconda uma 
implantação preferencial das amostras, 


— Dispositivo a 2 dimensões: 


Este princípio de composição directa de estimação elementar estende-se a configuraçõs a 2 dimensões: 


Seja, por exemplo, a estimação do teor médio de uma superficie 
poligonal S constituida por N áreas elementares s, reconhecidas cada 
uma por amostra central, 

Com raciocínio idêntico aos anteriores, obtêm-se: 


1 N e. 
Cy E = 18, 6 
Ss 
j=1 
Sendo s, = s, vem: 
=» Z 
“EN 


Seja ainda, por exemplo, a estimação do teor médio numa super- 
fície S constituida pela reunião de N rectângulos de área s = L 1, 
cada um deles reconhecido pela sua mediana L. Então: 


a 

q ec 1.7 

g i Cs 
Sendo 1 = 1, vem 

: 1 2 2 

é L, Ti 
À F | 
(ZL) i 
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— Dispositivo a 3 dimensões: 


Para configurações a 3 dimensões, o princi- 
pio aplica-se de modo idêntico: por exemplo, para 
um volume V constituido pela reunião de N volu- 
mes v, = 5 h, cada um deles reconhecido pela 
sua secção mediana, será: 


2 ! LX E 
a PR 
ve 1 
NS o = pps Eu s; o 


A aproximação, neste caso, consistiu em supôr os erros de estimação elementares de cada volume como 
independentes entre si. 


b) Composição de termos 


A informação disponível para o reconhecimento de um fenómeno a 3 dimensões está, frequentemente, 
concentrada sobre linhas preferenciais (sondagens, grupos de sanjas, galerias, chaminés, etc. numa jazida 
mineral), ou sobre planos (os diferentes níveis de implantação das galerias, por exemplo). 

A estimação global do teor médio de um volume (uma massa mineral) pode ser decomposta, por exem- 
plo, em três etapas, correspondendo aos três termos da variância de estimação: 


— Estimação das galerias pelo conjunto das amostras (rouços) disponíveis: é o termo de Inha. 

— Estimação de cada nível pelas galerias (supostas perfeitamente reconhecidas): é o termo de plano 
ou de secção. 

— Estimação de cada «fatia» da massa pelo seu plano médio (nível), também suposto perfeitamente 
reconhecido: é o termo de fatia ou de camada. 


Fatia 


Nivel 


Galeria 


a a a 


Parse Sigea 


O ————— = 


Termo de Fatia Termo de Plano Termo de linha 
2 q. 2 
Uf p 24 


Se se admitirem os três erros de estimação, correspondentes a estas três etapas, considerados como 
independentes, a variância de estimação global da massa é constituída pela soma de três termos: 


2 2 2 2 
p= mto To; 
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Para cada um dos termos, e atendendo ao que já se referiu, será: 


2, 1; (para N amostras totais regularmente espaçadas nas 
=>— q 
E de galerias, sendo o: a variância de estimação elementar 
do ponto central ao segmento e) 
1 M 
ci=—— E so? (o; é a variância de estimação do rectângulo s, pela 
cá j=1 4 sua mediana L, suposta perfeitamente reconhecida) 
a 1 as 2 E ss | 
Cf — EE ps ME ER (7,, é à variância de estimação do paralelipípedo v, pelo 
k=1 seu plano meridiano s,, suposto perfeitamente reconhe- 


cido) 


Esta composição de termos apresentada pode, evidentemente, ser aplicada de modo análogo a outras 
configurações. 


45-— O problema do suporte: regularização e acumulação 


Uma V.R. é sempre definida num certo suporte caracterizado pelo seu volume, forma e orientação. 
Portanto, munca podemos conhecer, na prática, os valores pontuais z(x,) das variáveis estudadas, mas sim 


valores regularizados por uma certa função de ponderação; por exemplo, o valor médio z, (x) no volume v 
amostrado: 


É 
ZQ(X) = E [ z(y) dy 


O problema prático que se põe é o de, a partir do variograma experimental regularizado (o único a 
que temos acesso), reconstituir o variograma pontual, visto ser este que deve entrar no cálculo das variâncias 
e covariâncias, 

Em certos casos, o suporte da regularização é bastante pequeno, podendo considerar-se pontual, e então 
é possível trabalhar directamente sobre o variograma regularizado. Quando tal não for possível, terá de se 
trabalhar sobre variáveis regularizadas. 

Define-se a regularizada de uma função z(x) por uma função de ponderação p(y), como o produto 
de convolução (*): 


+ 00 É 
Z (x) = [ z(x+y) p(y) dy =2*p 
a 
w 
em que: p(y) = p(-y) é a transposta de p(yv). 


O valor médio precedente z (x) corresponde à função de ponderação proporcional à indicatriz geomé- 
trica do volume v: 


1 
Pix) =— Es) 
v 
com: 
e 1 seyev 
k e 
dai l O seygév 


Pode-se interpretar 24x) como uma realização particular da F.A. 2! x), dita regularizada da F.A. 
a(x) pela função de ponderação p. 


a) Variograma regularizado: 


Se Z(x) é uma F. A. estacionária de 2.º ordem e admite a covariância C(h), a regularizada Z 51X) [- 
também F. A. estacionária de 2.º ordem e admite a covariância: 


(+*) O produto de convolução de duas funções f e g é definido por: 
f*g= [ f(y) g(x-y) dy 
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=". 


Ch) = E [Z,(x) «A (x+h)] = 
=E[ fZMx+y) «P(y) dy. f Z(x+h+w) .p(w) dw] = 
= SJ E[Z(+y) «Zixth+w)].p(y).p(w) dy dw 
Fazendo a mudança de variáveis: w=y+u, teremos: 
Co th) = f f C(h+u).p(y).p(y+u) dy du = 
= ] C(h+u) du. f p(vt+Hu).p(v) dy 


Mas este 2.º integral é, por definição, o covariograma transitivo (*) da função de ponderação p(v); 
logo: 


C(h)=f C(h+u).P(u) d=C*P=C*P 


—— 


P=(p*'p)=p*p=p*p=P 


Isto é, a covariância da regularizada Z, obtém-se regularizando a covariância inicial pelo covariograma 
transitivo da função de ponderação. 


O variograma regularizado deduz-se então de C,(h) pela relação clássica: 
Yo(h) — C(O) — Ch) 


Pode-se também deduzir directamente Yo lh) do variograma pontual y(h), substituindo na expressão 
anterior e obtendo: 


Yolh) = f ylx+h).P(x) dx— ( y(x). P(x) dx 
Ou, em notação simbólica: y =y * P—A 
Sendo a constante A determinada pela condiçãoy (0) = O. 


Na prática, a ponderação p corresponde geralmente a uma média sobre o volume v do suporte, pelo 
que o variograma regularizado se pode interpretar como a variância de estimação da amostra v(x) pela 
amostra v(x+h): 


1 
um) =— E[Z, (x+h) — 2, 009] : 


Pode-se então usar a fórmula, mais prática, do cálculo da variância de estimação de uma média por 
outra: 


TM) =Y(v va) — y(vv) 


em que Y( V,V,,) representa o valor médio do variograma pontual y(x-y) quando os dois pontos de apoio 
percorrem, independentemente um do outro, um o volum v e o outro o seu transladado de h(v,,). 


Note-se que, para distâncias h » v, vem: y(V,V,,) = y(h), donde: 


44h) =y(h) — Y(v,v) 


(+) Define-se o covariograma transitivo de uma função p(v) como: 


P(u) = f p(y) «p(y+Hu) dy. 
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Quer dizer, para distâncias h >» v, o comportamento pontual do variograma conserva-se, a menos de 
uma translação para baixo de valor: 


y(v,v) = Ev) 
O alcance é aumentado de: a=a + dim(v) 


Demonstra-se também (V. por exemplo F. MUGE et al. (15) ) que a constante de pepita C se regu- 
lariza na razão inversa do volume de regularização, isto é: 


Why a | E 
Uvsa | Ty = 


to) = O 
v 


b) Acumulação 


Designa-se por acumulação a regularização correspondente à média calculada sobre possança constante. 
Considere-se, por exemplo, o caso corrente na prática de sondagens verticais, sendo Z(y) a F. A. pon- 


A 8 tual com variograma y(h) e Z,(x) o teor médio no comprimento 1 
o (constante) de um troço de sondagem: 
à ul 
Z, (x) = — f 20) dy; x representa as coordenadas (X, X,) no 
k: | 
plano horizontal da sondagem, e y a coordenada vertical. 


O variograma acumulado y,(h) numa direcção horizontal 


dd 
qualquer h escreve-se: 
1 7 eu = 
Y,(h) = = E[Z,(x) — 2, (x+h)] = (1, 1) — sell, 1) 


O símbolo y (1, 1) designa o valor médio do variograma 
pontual y(h) entre os dois lados paralelos de comprimento 1 (os 
troços das sondagens). 


Se o variograma pontual y(h) é isótropo no espaço a 3 dimensões, o variograma horizontal acumu- 
lado y,(h) é também isótropo mo espaço a 2 dimensões e tem por expressão: 


2 2 
l paia l ps = 
nO) = 1-9) v(vRty?) dy —— f0-9) v() dy = 4(141)-70,D = a(h,D-a(o, 1) 


Na prática, desde que h > 1, pode-se assimilar o variograma acumulado y(h) ao variograma pontual 
y(h) transladado para baixo de um valor igual à diferença entre os patamares (no caso de esquemas de 
transição), isto é: Yth) = y(h) — F (1). 

Note-se que é necessário distinguir o variograma regularizado (acumulado) por acumulação sobre uma 
possança 1, do variograma regularizado por troços de comprimento 1 empilhados, embora o suporte da regu- 
larização seja o mesmo nos dois casos («carottes» de comprimento 1, por exemplo). No caso da acumulação, 
os dois segmentos 1 e 1, são paralelos, enquanto no caso da regularização eles estão alinhados. 

Finalmente, é de referir que a acumulação permite deduzir as funções auxiliares a (n+1) dimensões 
das respectivas expressões a n dimensões. 


q 1 
Por exemplo: seja o rectângulo ABCD, no qual se define a função auxiliar: x(L,1) = A fu, 1) dx, 
sendo a calculado sobre valores pontuais. | 
Se projectarmos a figura no eixo dos xx, realizamos uma 
acumulação sobre a possença 1 A expressão de Xp (LO), calculada 
agora sobre valores acumulados, é: 


a pl e qd dx É pE O Dai = 
Rg dig to DE S 26x D > Su “Dal = 
= y(L,1) — (0,1) 
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De modo análogo se definem as correspondências (por regularização R): 


R 
v(L) T+ aq(L, 1) — a(0,1) 


R 
x(L) Ee ==. x(L,1) — y(0,1) 


R 
F(L)”-—> F(L,1) — F(o,1) 


Para 3 dimensões, o processo de correspondência é similar. 


5 — APLICAÇÕES 
5.1. Aplicação à optimização de um reconhecimento 


Uma aplicação interessante do conceito de variância de estimação é na optimização de campanhas de 
reconhecimento de jazidas minerais: conhecidos os custos unitários de amostragem, determinar o plano de 
amostragem que conduz ao menor erro para dado custo global fixo, ou então o que implica a menor despesa 
total para dado erro previamente fixado. 


Um exemplo simples permite concretizar este problema: suponha-se que se pretende reconhecer um 
filão estreito sub-vertical, de configuração aproximada mente rectangular LXH, por galerias de nível amos- 
tradas descontinuamente por rouços traçados a intervalos regu- 

lares. 


E Existe uma verba global P= 20X 10º unidades monetárias 


alocada à campanha de reconhecimento, e computam-se em 
P, = 5000 unidades e p =500 unidades os custos do metro de gale- 
ria e da abertura do rouço, respectivamente. Sejam: L=1000m, 


H=300m. 
H p— | —s Por uma análise estrutural previamente feita na região, 
CO apta ss sabe-se que a V.R. em estudo (o teor, por exemplo) segue um 


esquema de tipo esférico isótropo, com alcance a=36 m., 
Pretende-se então optimizar o plano de amostragem, deter- 
minando o afastamento h entre galerias e 1 entre amostras que 
conduzem ao menor erro de estimação do filão. 
Admitindo o princípio de composição dos termos de linha 


E e e de plano, a variância de estimação total será dada por: 


E —>>—>—— LA 1 h dd 
Das a 2 2 
Cp — 6 Ca) = A) 


( 


1 

m 
em que of é a variância de estimação de um segmento 1 pela sua amostra central, a a variância de estimação 
do rectângulo 1.h pelo seu segmento mediano 1, m=H/h o número de galerias a traçar, e n=L/1 o número 


de amostras a colher por galeria, 
Por outro lado, o custo total dos trabalhos de reconhecimento será: 


H 1 
P=mLp RAR RE, da 1) 


condição esta que é necessário respeitar. 


Como se sabe do cálculo das variações, o mínimo de E atendendo ao constrangimento P fixo é dado 
pelo sistema: 


ê ) 
qa (op +A P)=O 
21 | 


(0, +h P)=O 


P 
% » (multiplicador de Lagrange) 
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Designando por r a relação P,/P, o sistema escreve-se: 


h o : 
= a Met) A 


h 


RE E 
+ + — 
HP Hadh 


3 1 r SB 
nu a tido 


AR e 
“a Pe l 1) 


Ou, eliminando ).: 


P) É 2 


L á Hg | 
Es h.o )= 


Em ambas as equações, o 1.º membro só depende deh e o 2.º só de 1, pelo que é possível achar a inter- 
cepção das duas curvas no plano [h,1] que fornece o par de valores óptimo procurado. 

O cálculo analítico é simples, dado que se conhece o modelo do variograma válido para a regionali- 
zação. Admitindo, como hipótese razoável, quel<aeh>a, as variâncias de estimação escrevem-se: 


a Es PD=C— 
Cc, = — | = 

| x 65) 6a 
És quit 1 EPMD ass des 
Gs — mi Bee ci = 

p= 2x 67) 4L 5º 6Lh TOL?h? 


) +€ 


Substituindo, valores e resolvendo o sistema, obtém-se o par óptimo: 


1=6,70m 
h=76,8m 


O que leva a traçar é galerias distanciadas de cerca de Tôm, e em cada uma delas abrir 150 rouços, 
distanciados de aproximadamente 6,67 m, ou seja, colher um total de 600 amostras. 

O custo total será de 20,3. 10º unidades monetárias (devido às aproximações feitas). Como era de 
esperar, a maior parte desta verba refere-se à traçagem das galerias de nível, 


5.2 — Aplicação ao cálculo de reservas 


Outra aplicação prática da estimação geoestatística é o cálculo da reserva mineral de jazidas. Este 
problema será retomado e desenvolvido no capítulo sobre krigagem, todavia apresentam-se desde já algumas 
noções relativas aos erros de estimação. 

Assim, numa estimação global podem-se distinguir dois tipos de erros: 


— um erro geométrico, ligado à ignorância dos limites exactos da mineralização; 
— um erro de qualidade, ligado à estimação dos teores e outras características da mineralização. 


O erro de qualidade afecta, por exemplo, a estimação dos teores médios, isto é, das quantidades de 
metal existentes. O erro geométrico, muito mais complexo de definir, afecta a estimação da geometria do 
jazigo e, portanto, as estimações de volumes e tonelagens de minério e também as quantidades de metal. 
a) Erro de qualidade: 


Considera-se, numa 1.º fase, que a geometria e volume V do campo são conhecidos. Pretende-se esti- 


1 

mar a característica média Z, = ae JS Z(x) dx no campo (por exemplo o teor médio), sendo Z(x) o teor 
v 

pontual, através de um estimador Z* (que, se for o estimador de krigagem, é óptimo), e que na maior parte 
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das vezes é uma combinação linear de estimadores parciais (de painéis, blocos, pisos, etc.). Supõe-se, eviden- 
temente, o campo estruturalmente homogêneo, de variograma y(h). 
Então, o estimador global será dado por: 


1 
e ma a E com V= 44 Y, 
a i 


A variância de estimação global será calculada do modo já atrás referido. 


b) Erro geométrico 


O erro geométrico pode estar directamente ligado a hipóteses de ordem geológica (génese da jazida, 
tectônica, etc.), e então tem um carácter determinista pouco favorável à sua quantificação, ou pode apresen- 
tar um carácter parcial ou puramente aleatório, que é então possível de ser traduzido num modelo quantifi- 
cável utilizando, por exemplo, a simulação, 

G. MATHERON elaborou um método puramente probabilístico, formalizado na Teoria Transitiva, que 
permite a dedução de expressões para a variância de estimação de espaços a uma, duas ou três dimensões. 

A V.R. neste caso é a variável em tudo ou nada, indicatriz geométrica K(x), definida por: 


1 se xXxedJ 
0 se xgéJ 


O jazigo J é definido pela sua massa mineralizada (total ou, por exemplo, acima de um dado teor de 


corte), de volume: 
+ co 
V- f Ki(x) dx 


o — 20 


Ki(x) = 


Em teoria transitiva define-se a função estrutural covariograma transitivo como: 


+ 50 
P(h) = f Ki(x,. K(x+h) dx 
“ —So 
Note-se que V=P(O), como é evidente. 
O erro de estimação geométrico é dado por: 


.) + 0 + DO 
opta) = a ss P(p.a) - | P(h) dh 
p=  —0 
em que a é o volume da malha (Ver demonstração em G. MATHERON (11)). 

Atendendo a que, quaisquer que sejam a forma e a grandeza de V,P(h) tem um comportamento linear 
na origem, o primeiro termo do desenvolvimento da expressão de cela), quando a tende para zero, pode 
ser calculado, o que permite obter fórmulas aproximadas de fácil manipulação. 

Assim, para 2 dimensões, vem para a variância relativa: 


o, 1 1 N, 
ie E — N, + 0,06 — [com n > 10] 
Sº nº 6 N, 


Sendo 2N, e 2N, (com N, < N,) os números dos elementos paralelos aos dois lados da malha (a, a,) 
que constituem o perímetro da reunião S* das n zonas de influência positivas. Note-se que o perímetro de 
eventuais lacunas no interior de S* entra no cálculo de N, e N., mas existindo partes do perímetro perfei- 
tamente reconhecidas, estas não são contadas (por exemplo, limites de falhas geológicas, etc.). 

Assim, sendo o estimador da reserva (em minério ou metal) medido pelo produto ou quociente de 
vários estimadores parciais (superfície, possança da camada, teor, densidade do minério, por exemplo), o 
erro global de estimação será calculado atendendo às regras de estimação de produtos e quocientes, que 
apresentam normalmente problemas dificeis de convergência. 

Para variáveis seguindo leis lognormais, demonstra-se que: 


2 2 2 
Se: Pt = A*,B*, então: Pp Ta TB TA Tg 
=="— + 2 j EF 
Pi A Bº AB A B 
A* 2 E 2 
Se: Qr=—, então: “o CA Sp " Ta 4 
* — TT — + —— — —— .- 
Q: A? Bº p AB A EB 
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E((A-A*) (B-B*)] ' eo 
-—— wo coeficiente de correlação dos dois erros de estimação. 


Sendo: 
RB CA: Op 

Estas expressões são, em geral, aplicáveis sem grande erro mesmo para variáveis não lognormais. 

Em L. CORTEZ et al. (1), apresenta-se um exemplo prático da estimação da reserva de um jazigo 
do tipo camada, comparando os erros obtidos por diversos métodos. 

Apresenta-se seguidamente outro exemplo prático de um cálculo global de reservas de um jazigo tam- 
bém do tipo camada, reconhecido por malha de amostragem regular. Este exemplo será bastante simplifi- 
cado, visto se pretender apenas esquematizar os passos fundamentais do cálculo e comparar a precisão das 
hipóteses de aproximação referidas atrás. 

Considere-se o jazigo estratiforme hematítico em cuja planta se assinalam as sondagens executadas 
para reconhecimento, com malha regular de 200m X 100 m. 


Sondagens positivas: x 
PLANTA DO JAZIGO Sondagens negativas: o 


Em cada sondagem central positiva, mediu-se a espessura mineralizada e analisou-se o teor em Fe. 
Os resultados indicam-se no Quadro seguinte: 


N.º Sondagem Espessura (m) Teor (% Fe) Acumulação (m %) 


1 0,8 51,2 40,96 
2 1,4 54,7 76,58 
3 2,0 55,4 110,80 
4 0,7 47,3 33,11 
5 1,3 48,9 63,57 
6 2.0 58,2 106,40 
7 1,8 54,0 97,20 
8 1,8 50,5 80,80 
9 1,5 54,0 81,00 

10 2,2 50,1 110,22 
11º 2,4 49,3 118,32 
12 0,9 46,0 41,40 
13 0,9 47,5 4275 
14 2.0 50,5 101,00 
15 1,6 49,9 79,84 
16 2,7 48,6 131,22 
17 0,8 47,0 37,60 
18 1,8 45,2 81,36 
19 1,6 49,5 79,20 
20 2,0 47,4 94,80 
21 1,0 45,3 45,80 
Valor médio 1,57 49,79 78,73 
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Suponha-se que se pretende calcular a tonelagem total de minério e a quantidade total de metal (ferro) 
contido no jazigo, bem como o seu teor médio. 


A tonelagem total de minério será dada por: 


+) 


T=B,o. 


em que S é a área total mineralizada, ; a densidade média do minério e h a espessura média mineralizada. 
O seu estimador será: 


T* =S*.p*.hº 
Vamos supor que a densidade média se pode estimar sem erro, o que aliás é hipótese vulgar na prá- 


tica. Seja: p* = 2,9ton m:3. 
Como estimador da superfície mineralizada, toma-se a união das 21 áreas de influência de cada son- 


dagem positiva. Isto é: 
S* = 21 x 200 x 100 = 420 X 103 m? 


Para uma estimação global, é em geral suficiente tomar como estimador da V.R. (neste caso, a espes- 
sura mineralizada) a média (ponderada pelas áreas de influência, no caso de malhas irregulares) dos valo- 


res amostrados. Ou seja, neste caso: h* = 1,57m. 
Donde: Tº* = 420 X 29 X 1,57 X 103 = 1.912.260 ton 


Este estimador está, evidentemente, afectado de um erro, que será medido pela sua variância relativa: 


o ps ” 


Dado que estas duas valláveis são manifestamente independentes e se admite que o erro na densidade 
é desprezável. Note-se que: 


é a variância do erro geométrico, 


é a variância do erro de qualidade. 


Usando a expressão da teoria transitiva já indicada, e atendendo a que n=21, N,=N,—9, vem para 
o erro geométrico cometido na estimação da superfície mineralizada a variância: 


2 
Em 1 qd N, 1/9  0,06x9 
a — o N, + 0,06 en) e ( RR a 

S2 n2 16 N 21º 16 9 


= 0,00463 


à 

O cálculo do erro de qualidade será feito através da variância de estimação, utilizando os princípios de 
aproximação já citados e, dado que a malha de reconhecimento é regular, recorrendo a ábacos. E óbvio que 
para tal é necessário traçar os variogramas experimentais da V. R. estimada e ajustar-lhe um modelo teórico. 

Dado que a pequena extensão do jazigo que se escolheu para este exemplo não permite o traçado de 
variogramas suficientemente representativos e coerentes, vamos supor que o jazigo está englobado numa 
mineralização mais vasta, cujo estudo estrutural revelou um variograma da espessura que é ajustado por 
um esquema esférico isótropo com amplitude a = 400m, Co = 0,20 e C = 0,85. 

A variância do erro cometido ao tomar a média das espessuras amostradas como estimador da espes- 
sura mineralizada pode ser avaliada de três modos: 


a) — pela variância de estimação dos rectângulos de influência através das suas amostras centrais; 


b) — pela decomposição em fatias com a direcção q, (de 100m de largura) reconhecidas através 
das suas linhas médias supostas conhecidas, sendo estas estimadas por extensão das amostras cen- 
trais aos seus segmentos de influência (de 200m); 
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poder 


c) — idêntico ao anterior, mas com fatias na direcção q, (de 200m de largura). 


Em qualquer dos casos, como já se referiu, é necessário admitir a independência do= erros para se 
fazer a sua composição. 
Far-se-à o cálculo pelos 3 processos, para efeitos de comparação de resultados: 


a) — A variância de estimação do rectângulo de 200m X 100m pelo seu ponto central é dada por: 


200 100 
s( mo) =02905C+€C 
400 400 


Donde ci=2X (0,22 0,85 + 0,20)- (0,295 x 0,85 + 0,20) = 0,323 
Supondo os erros de estimação dos 21 rectângulos independentes: 
>» 0,328 


p= —— 0,0154 


Ou seja, uma variância relativa: 


[o 
| — = 0,00463 + 0,00625 = 0,0109 


O que corresponde, ao nível de confiança de 95 %, a um erro relativo de: 


Sr 
2 


= 0,2086 = 20,86 % 


Ou seja, com 95 % de probabilidades, a tonelagem total está compreendida no intervalo: 
[1 531 363 ton, 2311157 ton]. 
b) — O erro de linhas é calculado pela variância de estimação aos segmentos de influência de 200m 
l 200 
| 


dos seus pontos centrais. Esta variância é dada no ábaco n.º 7: gp, para — = Oo. = 0,5 
; a 


= 0,13 C+C = 013 XxX 0,85 + 0,2 = 0,3105 
A variância de extensão dos 21 pontos centrais aos seus segmentos será evidentemente: 


0,3105 
21 


Fi 


o = = 0,0148 


O erro de fatias tem por expressão: 
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Visto as fatias terem comprimentos 1 diferentes. A variância de extensão da mediana ao rectângulo de 


10 
influência é dada no ábaco n.º 8, para L/a — oo = 0,25 e | variável consoante os comprimentos das linhas. 


Na tabela seguinte indicam-se os valores obtidos, estando as fatias numeradas de cima para baixo, 


Fatia Comprimento a 3 a 2 3 
Cá T l,/a Sm q SIT] + É oa 
1 800 2,0 0,00425 | 0,0036 2 304 
z 800 2,0 0,00425 0,0036 2 304 
3 600 1,5 0,0050 Ê 0,00425 1530 
4 1200 3,0 0,0035 0,0030 4 320 
3) 600 1,5 0,0050 0,00425 1530 
6 200 0,5 0,0125 0,0106 | 424 
z 4200 va = em 12412 
2 12 412 
Donde: o = -—-——— = 0,0007 
42007 


Então: q, =, +, = 0,0148 + 0,0007 = 0,0155 


Logo, a variância relativa vem; 


fia 0,0155 
DE — E 
ri 1,57º 


Para a tonelagem total, vem: 


Í 
Gg 
a — 0,00463 + 0,0063 = 0,01093 


Para o mesmo nivel de confiança de 95 %, vem o erro relativo de: 


Tr 
2 KT = 0,2091 = 20,91 % 


Note-se que a utilização deste processo de decomposição em termo de linhas e termo de fatias con- 
duz a um erro relativo quase idêntico ao processo de decomposição em rectângulos de influência. Note-se 
ainda que o efeito de pepita apenas entra no termo de linha, não entrando como factor aditivo no cálculo 
da variância de extensão do termo de fatia. 


100 
c) — Neste caso, para =>—— , o erro de linhas é: 


l 
a 


o. = 0,064 C+C = 0,2544 


5 0,2544 
o = ———— = 0,0121 
21 
dd L | 200 
O erro de fatias é obtido por uma tabela semelhante ao caso anterior, para —- =-—— | e sendo as 
a 400 


fatias numeradas da esquerda para a direita. 
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1 
2 
J 
4 
) 
6 
bh 
ui 


17.450 
comem 2 OUUÓS 


É 
Donde: q = — 
: 2100º 
Logo: cp = 0,0121 + 0,0040 = 0,0161 


= 0,00653 


aa e DE 0,0161 
seja. = — 
= 1,57? 


[2] 
Então: e — 0,00463 + 0,00653 = 0,0112 


Donde, ao nível de confiança de 95 %: 


É 


Como era de esperar, o erro obtido por composição de termos de linhas e de fatias é praticamente 
idêntico, independentemente da direcção em que se tomam as linhas. 

Para o cálculo do metal contido e do teor médio do jazigo, utilizar-se-ão as acumulações a = ht, 
onde t é o teor tal-qual, como V. R. 

Admite-se que o variograma das acumulações é bem ajustado por um esquema esférico isótropo, de 
amplitude idêntica (a=400m), C, = 0,10 e €C = 0,50. 

A quantidade de ferro total pode ser estimada por: 


E ot E 


Q =s.p .a = 420.10" x 2,0 XxX 0,788 = 959.784 ton Fe. O erro deste estimador terá por variância: 


Admite-se ainda que o erro na estimação da densidade média é desprezável, e que os dois erros 
(geométrico e de qualidade) são independentes. 

O valor de 04/82 já foi anteriormente calculado. A variância da acumulação média pode ser avaliada 
por qualquer dos processos atrás descritos. 

Utilizando o primeiro processo (que dá menor erro), obtém-se: 


o; =2X (0,22 X 0,50 + 0,10) — (0,295 X 0,50 + 0,10) = 0,1752 


» 01725 
d=—.=000%m 
| 
co 0,00821 
—— — = (0156 
= 0,725? 
2 
Ts 


Como: E = 0,00463 (já calculado), vem: 


TÉCNICA 451/452 135 


2 
[e] 
2 — 000463 + 0,0156 = 0,02023 


Q? 
Ou, ao nível de certeza de 95 %: 


Ta 
dE = 0,2845 — 28,45 % 


Ou seja, com 95 % de probabilidade, à quantidade de Fe está contido no intervalo: [686.725 ton, 
1 232 843 ton]. 


O teor médio pode ser estimado por: 


E á (9º a 


ft — E 


é hº 


= 0,5015 — 50,15 % Fe 


Como, em geral, não é possível fazer a hipótese da independência dos erros de acumulação e de espes- 
sura, a estimação do erro do cociente será: 


à ESTA pace RE sa 
pr a h* a a h 
Dado que: 
e [0] — 
+. = 0:0158 ". = 61480 
hº a: 
) ex 
E 0,0223 o =='0:1400 
ds ho 


E o coeficiente de correlação entre as acumulações e as espessuras é: 


o = 0,946 


ah 
Vem: 


Gg —- 
-*º = 0,0156 + 0,0223 — 2 X 0,946 x 0,125 + 0,149 = 0,00266 


[ef - 
Ou ainda: 2 * = 0,1032 = 10,32 %. 


Isto é, com 95 % de probabilidade, o teor médio estã contido no intervalo [44,97 % Fe, 55,33 % Fel. 

Note-se que, nos cálculos efectuados neste exemplo muito simples, se fizeram dois tipos de aproxima- 
ções que afectaram a precisão da estimação dos erros: 

a) —admitiu-se o princípio da independência dos erros, na sua composição; 

b) — utilizaram-se ábacos no cálculo das variâncias de estimação elementares, 

A comparação dos valores obtidos pelos diferentes processos permite avaliar a influência das aproxi- 


mações feitas e verificar assim, num caso prático, a validade dos princípios de aproximação usados na 
prática da estimação geoestatística. 
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Novo sistema de armazenagem de informações 


A Philips apresentou em Eindho- 
ven um novo processo de armaze- 
namento de informações que per- 
mite registar o conteúdo de 500 000 
páginas dactilografadas num disco 
plastificado de trinta centímetros. 


O conjunto apresenta-se sob a 
forma de um electrofone cujo dis- 
co, em vez de preto, é transparente 
e brilhante e onde o diamante foi 
substituído por um diodo laser 
cujo finíssimo feixe de luz forma, 
sobre a banda, uma mancha lumi- 
nosa com um milésimo de milime- 
tro de diâmetro. 


Embora não se verifique qual- 
quer contacto com o disco, É pos- 
sivel lê-lo estudando a luz que a sua 
superfície reflecte como um espe- 
lho. Com efeito, nas zonas de regis- 
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to da informacão — sob a forma de 
orifícios microscópios feitos na 
banda —, a luz não se reflecte tão 
bem e esta alternância de tons es- 
curos e brilhantes permite uma rá- 
pida descodificação, sendo o tempo 
de acesso a um da ordem de 
um quarto de segundo. Mas, para- 
lelamente com a forma de leitura, 
também a forma de registo elec- 
trónico assume especial importân- 
cia. No sistema lançado pela Phi- 
lips, basta aumentar a potência do 
diodo laser para fundir o material 
do disco e efectuar assim uma sé- 
rie de orifícios, característicos da 
informação a armazenar. O disco 
Philips tem capacidade para arma- 
zenar 10 mil milhões de «Bits», ou 
seja, o equivalente a uma enciclo- 
pédia em cem volumes, resultado 


que ultrapassa em grande medida 
os dos sistemas clássicos, ma opi- 
nião dos investigadores holandeses. 
A fórmula tem, porém, algumas 
desvantagens, pois, uma vez feita 
a gravação, não é possível apagá-la 
como acontece com os discos e ban- 
das magméticas. Os investigadores 
da empresa holandesa não desis- 
tiram, no entanto, de conseguir, 
até à comercialização deste produ- 
to, dentro de alguns amos, encon- 
trar entretanto solução para este 
problema, e lançar, no mercado, 
discos de todas as dimensões, mais 
aperfeiçoados ainda e portadores 
de dados, textos, comentários e 
imagens para todos os fins, 


[Extraído do «Boletim Mensal BEFNs, 
(12) Dezembro 1978]. 
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